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PRÓLOGO 


Este  texto  tuvo  su  origen  en  unos  apuntes  sobre  Ecuaciones  Diferenciales  para  los 
alumnos  de  la  Licenciatura  de  Matemáticas,  aunque,  a  lo  largo  de  estos  últimos  años, 
hemos  observado  que,  además,  resultaban  útiles  para  otras  carreras,  en  particular  para 
las  enseñanzas  de  Ingenierías  Técnicas  de  la  Universidad  de  La  Rioja.  Visto  que  estos 
apuntes  podían  ser  aprovechados  por  diversas  personas  con  diferentes  objetivos,  y  puesto 
que  podían  tener  un  público  no  demasiado  restringido,  nos  decidimos  a  darles  vida  en 
forma  de  libro. 

Los  métodos  clásicos  para  resolver  ecuaciones  diferenciales  son  importantes  pero 
difíciles  de  recordar.  Por  eso  nos  planteamos  escribir  algo  — en  principio,  los  apuntes 
antes  mencionados —  dedicado  a  ellos  con  exclusividad,  donde  se  pudiesen  encontrar  los 
métodos  fácilmente.  De  aquí  que  este  libro  110  contiene  nada  de  muchos  de  los  aspectos 
fundamentales  de  la  teoría  de  ecuaciones  diferenciales:  existencia  y  unicidad  de  soluciones, 
sistemas  de  ecuaciones,  integración  por  desarrollos  en  serie,  estabilidad,  . . . ,  por  citar  sólo 
unos  pocos.  Es  claro  que,  matemáticamente  hablando,  no  puede  plantearse  un  estudio  serio 
de  las  ecuaciones  diferenciales  sin  abordar  esos  temas,  pero  no  es  éste  el  objetivo  del  libro. 
Los  temas  que  aquí  se  tratan  pueden  explicarse  a  estudiantes  de  diversas  carreras  tal  como 
aparecen  desarrollados.  En  cambio,  el  estudio  de  la  existencia  y  unicidad  de  soluciones, 
por  ejemplo,  requiere  necesariamente  un  tratamiento  distinto,  ya  sea  más  práctico  o  más 
teórico,  dependiendo  del  tipo  de  personas  al  que  esté  destinado. 

El  libro  consta  fundamentalmente  de  tres  partes,  de  acuerdo  a  una  primera  clasificación 
general  de  la  ecuaciones  que  se  estudian:  ecuaciones  explícitas  de  primer  orden,  ecuaciones 
en  las  que  la  derivada  aparece  implícitamente,  y  ecuaciones  en  las  que  se  puede  reducir  el 
orden.  Cada  una  de  estas  partes  abarca  diversos  tipos  de  ecuaciones,  que  aparecen  en  lo  que 
liemos  denominado  “Apartados”,  y  que  hemos  numerado  consecutivamente  desde  1  hasta 
13.  Entre  estos  números  aparecen  a  veces  algunos  denotados  con  “prima”,  como  4'.  Alguien 
malintencionado  podía  pensar  que  tan  extraña  notación  respondía  simplemente  a  dejadez 
del  autor,  para  no  tener  que  renumerar  los  apartados  tras  haber  redactado  el  libro  en 
desorden.  No  es  éste  el  caso  (al  menos  en  estas  notas).  El  uso  de  “primas”  es  intencionado, 
y  quiere  significar  que  un  tipo  se  reduce  al  anterior  mediante  algún  mecanismo  en  forma  de 
cambio  de  variable.  Por  otra  parte,  todos  los  métodos  de  resolución  se  basan,  en  esencia, 
en  aplicar  transformaciones  diversas  hasta  llegar  a  una  ecuación  de  variables  separadas, 
cuya  resolución  requiere  sólo  calcular  integrales.  Así  pues,  no  tenía  sentido  utilizar  la 
denominación  1'  (o  sucesivas)  para  algún  tipo  concreto  de  ecuación,  puesto  que  lo  mismo 
podía  haberse  aplicado  a  la  mayoría.  Varios  de  los  tipos  que  se  estudian  se  subdividen  a 
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su  vez  en  subtipos.  En  todo  caso,  siempre  se  analizan  los  procesos  que  hay  que  seguir  para 
llegar  a  la  resolución,  a  veces  por  diferentes  caminos. 

Un  resumen  de  los  métodos  que  se  emplean,  para  recordarlos  de  un  vistazo,  es  lo 
que  aparece  en  lo  que  hemos  denominado  “Recetas”.  Estos  esquemas  permiten  clasificar 
fácilmente  las  ecuaciones  estudiadas  y  tener  una  rápida  indicación  de  cómo  abordar 
su  resolución.  Así  mismo,  con  cada  tipo  de  ecuaciones  se  muestra  un  ejemplo  típico 
completamente  resuelto. 

En  el  libro  aparece  una  pequeña  bibliografía  con  libros  exclusivamente  en  castellano.  Al 
contrario  que  en  muchos  otros  temas  de  matemáticas,  existen,  en  nuestro  idioma,  bastantes 
textos  dedicados  a  las  ecuaciones  diferenciales,  así  que  sólo  hemos  incluido  unos  pocos.  (La 
abundancia  de  libros  en  castellano  sobre  ecuaciones  diferenciales  se  debe,  en  opinión  del 
autor,  al  interés  del  tema  en  disciplinas  no  estrictamente  matemáticas.  Realmente,  en  los 
temas  más  puntuales  y  de  investigación,  esta  abundancia  ya  no  puede  considerarse  cierta.) 
Entre  las  obras  citadas,  no  hemos  considerado  necesario  indicar  cuáles  son  teóricas  y  cuáles 
se  dedican  fundamentalmente  a  la  resolución  de  problemas,  ya  que  nos  ha  parecido  que  sus 
títulos  son  bastante  descriptivos. 

Acaba  el  libro  con  un  apéndice  dedicado  a  los  métodos  de  resolución  de  integrales 
inmediatas  o  cálculo  de  primitivas.  Tal  como  ya  hemos  mencionado  anteriormente,  todas 
la  ecuaciones  que  aquí  estudiamos  se  intentan  reducir  a  ecuaciones  en  variables  separadas 
cuya  solución  se  expresa  por  medio  de  integrales.  Así  pues,  tal  recordatorio  puede  resultar 
claramente  de  interés  en  el  tema  que  estamos  tratando. 

Queremos  dejar  constancia  de  que  los  nombres  que  aparecen  en  el  índice  no  se 
corresponden  exactamente  con  los  títulos  que  hemos  ido  dando  a  los  diferentes  apartados. 
La  no  coincidencia  no  se  debe  a  descuido,  sino  que  ha  sido  pensada  conscientemente  para 
que,  cuando  alguien  se  encuentra  ante  una  ecuación  que  debe  resolver,  el  índice  le  permita 
una  rápida  identificación  del  tipo  que  se  trata,  y  dónde  se  puede  localizar  dentro  del  texto. 

Deseamos  así  mismo  justificar  la  falta  de  un  índice  terminológico  o  tabla  de  contenidos, 
que  quizás  alguien  pueda  echar  en  falta.  La  ventaja  que  tienen  tales  tipos  de  índices  es 
que  permiten  buscar  palabras  clave  clasificadas  alfabéticamente,  al  contrario  que  en  un 
índice  general  en  el  que,  obviamente,  los  apartados  aparecen  consecutivamente  según  el 
orden  en  el  que  se  abordan  dentro  del  libro,  y  en  el  que  muchos  términos  suficientemente 
descriptivos  pueden  no  estar  reflejados  o  ser  difíciles  de  localizar.  Es  opinión  del  autor  que 
casi  cualquier  libro  de  estudio  o  consulta  debería  llevar  un  índice  de  nombres,  así  que  no 
podemos  resistirnos  a  explicar  su  ausencia. 

Hay  que  tener  presente  que  éste  es  un  libro  pequeño  en  extensión,  dedicado  a  un  tema 
bastante  puntual,  con  un  índice  detallado,  y  cuyo  propósito  es  permitir  que,  cuando  nos 
encontramos  ante  una  ecuación  diferencial,  podamos  fácilmente  distinguir  su  tipo  para 
proceder  a  resolverla.  Así  pues,  no  parecía  demasiado  importante  algo  parecido  a  un  índice 
de  nombres,  ya  que  lo  que  interesa  al  lector  es  saber  identificar  el  tipo  de  una  ecuación  a  la 
vista  de  su  aspecto,  no  de  su  nombre,  que  es  fácil  que  quien  consulta  el  libro  no  conozca. 
Queremos  también  mencionar  la  dificultad  de  elaborar  un  índice  de  nombre  suficientemente 
completo;  esto  es  así  puesto  que,  aunque  muchos  de  los  tipos  de  ecuaciones  que  aquí  se 
estudian  sí  que  tienen  un  nombre  que  los  describe,  esto  no  es  así  en  todos  los  casos,  sino 
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que  muchas  veces  las  catalogamos  únicamente  por  su  aspecto.  Por  esta  razón,  además, 
muchos  de  los  títulos  de  los  apartados  son  meramente  descriptivos,  clasificando  el  tipo  de 
ecuación  mediante  una  fórmula.  De  todas  formas,  si  alguien  desea  buscar  una  ecuación 
por  su  nombre,  no  es  complicado  localizarla  en  el  índice  ya  que  éste  es,  necesariamente, 
pequeño. 

Tampoco  se  ha  incluido  un  índice  de  “recetas” ,  pues  siempre  aparecen,  como  mucho, 
un  par  de  páginas  después  de  cada  tipo,  luego  resultan  fáciles  de  localizar  a  través  del 
índice.  Lo  mismo  puede  decirse  de  los  ejercicios,  que  invariablemente  están  colocados  tras 
la  explicación  teórica  del  método. 

Aunque  el  libro  ha  sido  suficientemente  repasado,  y  ha  sido  ya  utilizado  como  apuntes 
fotocopiados  durante  varios  años,  la  experiencia  nos  muestra  la  práctica  imposibilidad  de 
evitar  que  se  deslice  alguna  errata.  En  este  aspecto,  es  de  destacar  que  todas  ellas  son 
debidas  al  autor  y  no  a  ningún  proceso  posterior  en  imprenta,  puesto  que  el  libro  ha  sido 
editado  directamente  a  partir  de  las  páginas  ya  impresas  suministradas  por  el  autor.  En 
su  confección  se  ha  utilizado  TjrjX,  a  cuyo  creador,  Donald  Knuth,  deseo  hacer  constar  mi 
gratitud  por  permitir  a  la  comunidad  matemática  (y  científica  en  general)  la  utilización 
de  tan  potente  y  útil  herramienta  destinada  a  elaborar  textos  de  gran  calidad  tipográfica. 
Lástima  que,  aún,  no  esté  lo  suficientemente  adaptado  para  escribir  en  lengua  no  inglesa. 

Así  mismo,  quiero  agradecer  a  mis  compañeros  del  Departamento  de  Matemáticas 
y  Computación  de  la  Universidad  de  La  Rio  ja  sus  sugerencias  y  correcciones  sobre  las 
versiones  preliminares  de  este  libro.  En  particular,  a  José  Luis  Ansorena,  José  Manuel 
Gutiérrez  y  Víctor  Lanchares,  cuyas  críticas  han  permitido,  sin  duda,  mejorar  el  texto. 
También  mi  reconocimiento  a  José  Javier  Guadalupe,  de  quien  aprendí  mis  primeras 
nociones  sobre  ecuaciones  diferenciales  hace  años,  cuando  estudiaba  en  el  entonces  Colegio 
Universitario  de  La  Rioja,  semilla  de  nuestra  actual  Universidad;  de  sus  apuntes  dictados 
en  clase  surgieron  parte  de  estas  notas,  que  se  han  ido  completando  durante  varios  años. 
Por  último,  a  mi  mujer,  María  José  Ramírez,  que  ha  soportado  mi  ausencia  durante  las 
múltiples  horas  que  he  dedicado  a  escribir  este  libro;  como  ahora  — sábado  a  las  ocho  de 
la  mañana — ,  que  duerme  en  la  habitación  de  al  lado  mientras  yo  doy  los  últimos  (¡ojalá!) 
retoques  al  texto. 


Juan  Luis  Varona 
Dpto.  de  Matemáticas  y  Computación 
Universidad  de  La  Rioja 
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Logroño,  febrero  de  1996 
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GENERALIDADES 


Desde  los  primeros  pasos  en  el  cálculo  diferencial,  de  todos  es  conocido  que,  dada  una 
función  y  =  f(x),  su  derivada  jé  =  f'(x)  es  también  una  función  que  se  puede  encontrar 
mediante  ciertas  reglas.  Por  ejemplo,  si  y  =  e-*3,  entonces  j)  =  —  3x2e~x3  o,  lo  que  es 
lo  mismo,  -p  =  —3 x2y.  El  problema  al  que  nos  enfrentamos  ahora  no  es  el  de  calcular 
derivadas  de  funciones;  más  bien,  el  problema  consiste  en:  si  se  da  una  ecuación  como 
^  =  —3 x2y,  hallar  de  alguna  manera  una  función  y  =  f(x)  que  satisfaga  dicha  ecuación. 
En  una  palabra,  se  desea  resolver  ecuaciones  diferenciales. 

La  forma  de  ecuación  diferencial  más  sencilla  que  puede  pensarse  es  -tL  =  f(x). 
Resolverla  consiste  en  encontrar  una  función  cuya  derivada  sea  f(x),  es  decir,  encontrar 
las  primitivas  (integrales  indefinidas)  de  f(x).  Por  tanto,  podemos  decir  que  los  métodos 
de  resolución  de  ecuaciones  diferenciales  constituyen  una  generalización  del  cálculo  de 
primitivas. 

Definición  1.  Llamamos  ecuación  diferencial  (E.  D.)  a  una  ecuación  que  relaciona 
una  función  (o  variable  dependiente),  su  variable  o  variables  (variables  independientes) ,  y 
sus  derivadas.  Si  la  ecuación  contiene  derivadas  respecto  a  una  sola  variable  independiente 
entonces  se  dice  que  es  una  ecuación  diferencial  ordinaria  (E.  D.  O.);  y  si  contiene  las 
derivadas  parciales  respecto  a  dos  o  más  variables  independientes  se  llama  ecuación  en 
derivadas  parciales  (E.  D.  P.). 


Ejemplos  de  ecuaciones  diferenciales  ordinarias  son 


dy 

dx 


-Ay 


2, 


(x  +  2 y)  dx  —  3y  dy  —  0 


y 


d2y 

dx2 


0; 


mientras  que 


du  du 

x^~  +y^~ =  u 

dx  dy 


y 


d3u 

dx3 

son  ecuaciones  en  derivadas  parciales. 


d2u  du 

dt2  dt 


(1) 

(2) 

(3) 

(4) 
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Otro  tipo  de  ecuaciones  que  pueden  estudiarse  son  las  ecuaciones  diferenciales  de 
retraso  (o  retardo ),  como  es  el  caso  de 

u'(t)  =  7  —  2u(t  —  3). 

Están  caracterizadas  por  la  presencia  de  un  desplazamiento  t  —  ¿o  en  el  argumento  de  la 
función  incógnita  u(t).  En  general,  son  más  difíciles  de  manejar  que  las  E.  D.  sin  retraso. 
No  nos  ocuparemos  aquí  de  ellas. 

Definición  2.  Se  llama  orden  de  la  ecuación  diferencial  al  orden  de  la  derivada  o  derivada 
parcial  más  alta  que  aparece  en  la  ecuación. 

Así,  por  ejemplo,  las  ecuaciones  (1)  y  (3)  son  de  orden  1,  (2)  es  de  orden  2  y  (4)  de 
orden  3. 

En  lo  que  sigue  nos  preocuparemos  sólo  de  ecuaciones  diferenciales  ordinarias  y,  como 
no  habrá  lugar  a  confusión,  las  denominaremos  simplemente  E.  D.  Por  lo  general,  salvo 
que  el  contexto  nos  indique  otra  notación  (o  ésta  provenga  de  los  cambios  de  variable  que 
efectuemos),  utilizaremos  x  para  denotar  la  variable  independiente  e  y  para  la  variable 
dependiente. 

Definición  3.  Decimos  que  una  ecuación  diferencial  (de  orden  n )  está  expresada  en  forma 
implícita  cuando  tiene  la  forma 


F(x,y,y',...,y^)  =  0 

siendo  F  una  función  F:fl  C  Rn+2  — >  R  con  O  un  subconjunto  (generalmente  abierto) 
de  Rn+2.  Y  decimos  que  está  expresada  en  forma  explícita  cuando  tenemos 

V{n)  =  f(x,y,y’,...,y^n~l)) 

con  f:D  C  Rn+1  — >  R  una  función  deñnida  en  un  subconjunto  D  (generalmente  abierto) 
de  Rn+1. 


Una  clase  importante  de  E.  D.,  bien  estudiada  y  con  buenas  propiedades,  es  la 
siguiente: 


Definición  4.  Se  dice  que  una  ecuación  diferencial  es  lineal  si  tiene  la  forma 


dny 

an^x)¿~d  +  an-i(x ) 


dn~ly 

(te™-1 


H - f  ai  (a )^|  +  a0(x)y  =  g(x); 


y  se  llama  lineal  homogénea  si,  además,  g(x)  =  0.  Dada  una  ecuación  lineal,  su 
correspondiente  ecuación  lineal  homogénea  en  la  que  se  ha  hecho  g(x)  —  0  se  denomina 
lineal  homogénea  asociada.  Una  ecuación  que  no  es  lineal  se  dice  no  lineal. 


Nuestro  objetivo  es  resolver  ecuaciones  diferenciales,  esto  es,  encontrar  sus  soluciones. 


Generalidades 
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Definición  5.  Decimos  que  una  función  y  —  ip(x)  definida  en  un  intervalo  I  (es  decir, 
tp:I  el  — »  M )  es  solución  de  una  ecuación  diferencial  en  el  intervalo  si,  sustituida  en 
dicha  ecuación,  la  reduce  a  una  identidad.  (En  otras  palabras,  si  satisface  la  E.  D.)  Una 
E.  D.  se  dice  resoluble  (o  integrable)  por  cuadraturas  si  su  solución  es  expresable 
mediante  integrales. 

En  general,  la  solución  de  una  ecuación  diferencial  de  orden  n  dependerá  de  n 
parámetros.  Pero  incluso  de  esta  forma  pueden  no  obtenerse  todas  las  soluciones  de  una 
E.  D.  Por  ejemplo,  cuando  tenemos  una  familia  uniparamétrica  de  soluciones  de  una  E.  D., 
una  sencilla  interpretación  geométrica  nos  muestra  que  también  la  envolvente  de  la  familia 
de  curvas  (si  existe)  es  solución  de  la  E.  D. 

A  continuación,  nos  dedicaremos  a  explicar  los  diversos  métodos  clásicos  de  resolución 
de  E.  D.  No  liaremos  hincapié  en  el  intervalo  de  definición  de  las  soluciones,  ni  efectuaremos 
un  estudio  detallado  de  la  rigurosidad  de  los  métodos  empleados  que,  en  esencia,  descansan 
siempre  en  la  regla  de  la  cadena  y  los  teoremas  de  la  función  inversa  y  de  la  función 
implícita.  No  nos  detendremos  nunca  en  comprobar  las  hipótesis  de  estos  teoremas,  sino 
que  supondremos  en  todo  momento  que  las  funciones  que  aparecen  en  los  métodos  descritos 
son  lo  suficientemente  “buenas”,  o  están  lo  suficientemente  restringidas  en  su  dominio,  para 
que  siempre  se  satisfagan  las  hipótesis  necesarias.  Tampoco  nos  preocuparemos  en  exceso 
de  saber  si  hemos  obtenido  todas  las  soluciones;  a  este  respecto,  en  algunos  casos  nos 
interesaremos  por  las  soluciones  singulares  de  una  E.  D.,  como  puede  ser  la  envolvente  de 
una  familia  de  soluciones. 

Nos  apresuramos  a  señalar  que  las  fórmulas  generales  que  aparecen  como  solución 
de  diversos  tipos  de  ecuaciones  no  deben  memorizarse;  más  bien,  el  procedimiento  debe 
desarrollarse  completo  cada  vez.  Para  ello  bastará  recordar  unos  cuantos  puntos  esenciales 
que  destacamos  en  las  cajas  de  texto  que  hemos  denominado  “recetas”. 

Por  último,  comentar  que,  en  los  ejemplos  que  nos  aparecerán,  el  lector  puede 
entretenerse  en  representar  gráficamente  las  soluciones  de  las  E.  D.  planteadas,  al  menos 
en  los  casos  más  sencillos  o  efectuando  un  simple  bosquejo  de  su  apariencia.  No  pensemos 
en  esto  como  una  pérdida  de  tiempo,  pues  ayuda  a  comprender  la  naturaleza  del  problema 
y  de  sus  soluciones. 


ECUACIONES  EXPLÍCITAS  DE 
PRIMER  ORDEN 


Son  las  que  tienen  la  forma 


y'  =  f(x,y). 


APARTADO  1. 

Variables  separadas. 

Si  tenemos  la  E.  D. 

g{x)  =  h(y)y\ 

formalmente,  podemos  poner  g(x)  dx  —  h(y)  dy\  si  suponemos  que  G  es  una  primitiva  de 
g  y  H  una  de  h,  tendremos  G'(x)  dx  —  H'(y)  dy  e,  integrando,  G(x)  —  H (y)  +  C,  que  es 
la  solución  general  de  la  ecuación. 

Expliquemos  con  un  poco  más  de  rigor  por  qué  funciona  el  método:  Sea  y  =  <p(x)  una 
solución  de  la  E.  D.,  es  decir,  p(x)  debe  cumplir  g{x)  =  h((p(x))<p'  (x) .  Pero  como  H  es 
una  primitiva  de  h,  por  la  regla  de  la  cadena  se  tiene  g(x)  —  h(p(x))p'  (x)  =  ( H  o  p)'(x). 
Integrando,  G(x)  =  (Ho<p)(x)  +  C  (lo  que  antes  hemos  expresado  como  G(x)  =  H(y)  +  C), 
de  donde  p(x)  =  H~1(G(x)  —  C ). 

En  los  pasos  anteriores,  está  justificado  emplear  la  regla  de  la  cadena  cuando  p  y  H 
son  derivables,  lo  cual  es  cierto  sin  más  que  suponer  que  h  sea  continua.  Y  finalmente,  para 
poder  despejar  ip  mediante  el  uso  de  II  1  bastaría  con  exigir  además  que  h  no  se  anulara 
en  el  intervalo  de  definición,  con  lo  cual,  como  H'  —  h  ^  0,  H  es  creciente  o  decreciente 
luego  existe  II ~ 1  (en  otras  palabras,  como  la  derivada  de  H  no  se  anula,  el  teorema  de  la 
función  inversa  nos  asegura  que  existe  H~1). 

Las  ecuaciones  en  variables  separadas  son  las  más  sencillas  de  integrar  y,  a  la  vez, 
las  más  importantes,  ya  que  cualquier  otro  método  de  resolución  se  basa  esencialmente  en 
aplicar  diversos  trucos  para  llegar  a  una  ecuación  en  variables  separadas.  En  ellas  hemos 
visto,  con  todo  rigor,  qué  hipótesis  hay  que  imponer  para  que  el  método  que  conduce 
a  la  solución  esté  correctamente  empleado,  y  cómo  se  justifica  el  funcionamiento  del 
proceso.  A  partir  de  ahora  no  incidiremos  más  en  estos  detalles  que,  aunque  importantes, 
sobrecargarían  la  explicación.  El  lector  puede  detenerse  mentalmente  a  pensar  en  ellos, 
justificando  adecuadamente  los  pasos  que  se  efectúen. 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


En  cualquier  caso,  conviene  recordar  que  la  expresión  ^  es  simplemente  una  útil 
notación  para  designar  la  derivada  de  y  respecto  de  x.  no  un  cociente  de  dy  dividido 
por  dx\  ni  dy  ni  dx  tienen  entidad  en  sí  mismas.  Esta  notación  se  emplea,  no  porque 
sí  ni  para  introducir  confusión,  sino  que,  al  contrario,  se  usa  porque  es  consecuente  con 
los  enunciados  de  varios  importantes  resultados.  Ya  liemos  visto  cómo  resulta  adecuada 
a  la  hora  de  recordar  cómo  resolver  ecuaciones  en  variables  separadas  g(x')  =  h(y  ') 
descomponiendo  g(x )  dx  —  h(x)  dy  (como  si  ^  fuese  realmente  una  fracción)  e  integrando 
ambos  lados  de  la  expresión  anterior.  Pero  no  sólo  aquí  se  pone  de  manifiesto  la  utilidad 
de  esta  notación.  Por  ejemplo,  el  teorema  de  la  función  inversa  prueba  (con  las  hipótesis 
adecuadas)  que,  cuando  y  es  una  función  de  x,  si  se  despeja  x  como  función  de  y  se  cumple 

dx  1  1 

dy  dy/dx  y'{x)' 

es  decir,  se  produce  un  comportamiento  similar  a  si  estuviéramos  operando  con  fracciones. 
Análogamente,  si  tenemos  que  z  es  una  función  de  y  y,  a  su  vez,  y  una  función  de  x.  la 
regla  de  la  cadena  establece  que  la  derivada  de  la  función  compuesta  z(x)  es 

dz  dz  dy 
dx  dy  dx ’ 

que  es  como  si  simplificáramos  dy  en  los  supuestos  cocientes  de  la  derecha.  Esto  permite 
usar  las  notaciones  del  tipo  ^  y  su  comportamiento  como  si  fuesen  fracciones  como  regla 
nemotécnica  de  los  resultados  anteriores. 


RECETA  1.  Variables  separadas. 

Son  de  la  forma 

g(x)  =  h(y)y'. 

Formalmente,  se  separa  g(x)  =  h(y)  jé  en  g(x)  dx  =  h(y )  dy  y  se 
integra. 


Ejemplo  1.  Resolver 


+  (sen  x)y  —  0. 
dx 


Despejando,  ^  =  —(sena;)  dx  e,  integrando,  loga/  =  cosa;  +  C.  es  decir,  y  =  ecosx+c , 
Sin  más  que  tomar  K  =  ec  encontramos  las  soluciones  y  —  Kecosx.  Fijarse  que,  en 
principio,  parece  que  K  tiene  que  ser  positiva;  pero  en  realidad  la  integral  de  —  es  log  \y\, 
lo  que  nos  llevaría  a  soluciones  con  valores  negativos  de  K.  Por  último,  notar  y  =  0  (es 
decir,  tomar  K  —  0)  también  es  claramente  una  solución  de  la  E.  D.,  aunque  no  se  obtiene 
con  el  método  seguido.  Así  pues,  la  solución  general  de  la  E.  D.  es  de  la  forma  y  —  Kecosx 
con  K  G  M. 


E.  D.  explícitas  de  primer  orden 
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APARTADO  2. 

Ecuación  de  la  forma  y'  —  f(ax  +  by ). 

Si  a  =  0  o  b  =  0,  la  ecuación  es  separable.  En  otro  caso,  efectuemos  el  cambio  de 
función  y(x)  por  z(x)  dado  por  z  —  ax  +  by,  de  donde  z'  —  a  +  by'  y,  por  tanto,  y'  =  ^-F^. 
Entonces,  sustituyendo  en  la  E.  D.  obtenemos  =  f(z ),  es  decir,  z'  —  a  +  bf(z ),  que 
es  de  variables  separadas.  La  escribimos  como 


dz 

X  a  +  bf(z)  ’ 

con  lo  que,  integrando,  x  =  f  (a  +  bf(z))~1dz  —  <f>(z,C).  Así  pues,  las  soluciones  de  la 
E.  D.  de  partida  serán 

x  =  4>{ax  +  by,  C ), 

de  modo  que  hemos  encontrado  y  como  función  de  x  expresada  en  forma  implícita. 


RECETA  2.  Ecuación  de  la  forma  y'  —  f(ax  +  by). 

El  cambio  de  función  y(x)  por  z(x)  dado  por  z  =  ax  +  by  la 
transforma  en  una  de  variables  separadas. 


Ejemplo  2.  Resolver 

y'  -  exey  =  -1. 

Tenemos  y'  +  1  =  ex+y ,  con  lo  que  si  efectuamos  el  cambio  de  función  dado  por  la 
sustitución  z  —  x  +  y,  la  ecuación  queda  transformada  en  z'  —  ez,  es  decir,  dx  =  e~z  dz, 
ecuación  en  variables  separadas  cuya  solución  es  x  —  —e~z  +  C.  Volviendo  a  las  variables 
iniciales,  C  —  x  =  e~x~y ,  de  donde  log(C  —  x)  =  —  x  —  y,  y  por  tanto  la  solución  de  la 
E.  D.  de  partida  es  y  —  —  log(C  —  x)  —  x.  (Observar  que  no  nos  liemos  preocupado  — ni  lo 
haremos  de  aquí  en  adelante —  de  poner  módulos  cuando  al  calcular  una  integral  aparece 
un  logaritmo.  El  lector  podría  analizar  estos  casos  con  mucho  más  cuidado.) 


APARTADO  3. 

Homogéneas. 

Supongamos  que  tenemos  la  ecuación 


Métodos  clásicos  de  resolución  de  E.  D.  O. 


Para  resolverla,  hacemos  el  cambio  de  función  y(x)  por  u(x)  mediante  u  —  Así,  derivando 
y  —  ux  tenemos  y'  =  u'x  +  u,  es  decir,  u'x  +  u  =  f(u).  Esta  ecuación,  que  podemos  poner 
como  u'x  =  f(u )  —  u,  es  de  variables  separadas.  Vamos  a  solucionarla: 

•  Si  f(u)  ^  u,  podemos  escribir  f  e,  integrando,  ./'  =  log(§). 

Despejando  x  obtenemos  x  —  Ce^u^  con  (¡){u)  —  f  f^\_u  ■  Por  tanto,  las  curvas  con 
ecuaciones  paramétricas 

í  x  =  Ce+W 
\y  =  Cue 

son  solución  de  la  ecuación  diferencial  para  cada  C  e  R.  (Esto  constituye  una  familia 
de  curvas  liomotéticas:  una  curva  se  obtiene  de  otra  mediante  una  homotecia,  es  decir, 
multiplicando  los  valores  de  x  e  y  por  una  constante.)  A  veces,  es  conveniente  expresar 
estas  soluciones  de  otras  formas.  Siempre  puede  ponerse  x  =  Ce^y/X\  solución  dada 
mediante  una  función  implícita.  Y,  cuando  en  x  —  Ce*™  se  logra  despejar  de  alguna 
forma  u  =  H(x,C),  la  solución  de  la  E.  D.  queda  mucho  más  sencilla:  y  —  xH(x,  C ). 

•  Supongamos  ahora  que  existe  algún  uq  tal  que  f(uo)  =  Uq.  En  este  caso,  es  inmediato 
comprobar  que  la  recta  y  =  UqX  es  solución:  y'  —  uo  =  f(uo)  =  /(|),  luego  se  satisface 
la  ecuación  diferencial.  Este  tipo  de  soluciones  que  no  se  obtienen  con  el  procedimiento 
general  suelen  denominarse  soluciones  singulares. 

Nota:  En  general,  una  función  h(x,y )  se  dice  homogénea  de  grado  a  si  h(Xx,Xy)  = 
Xah(x,y).  Es  inmediato  comprobar  que  una  E.  D.  de  la  forma 

P(x,  y)  dx  +  Q(x.  y)  dy  —  0 

con  P(x,y )  y  Q(x,y)  funciones  homogéneas  del  mismo  grado  es,  efectivamente,  una 
ecuación  diferencial  homogénea  (despejar  y'  —  Y  extraer 

A  =  x  de  P  y  Q).  De  aquí  proviene  el  nombre  de  este  tipo  de  ecuaciones. 


RECETA  3.  Homogéneas. 

Son  de  la  forma 


Se  hace  el  cambio  de  función  y(x)  por  u(x)  mediante  y 
transformándose  así  la  E.  D.  en  una  de  variables  separadas. 


ux, 


,  2  xy  -  y 2 

y  = - 2 — 1 

xz 


Ejemplo  3.  Resolver 


E.  D.  explícitas  de  primer  orden 
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Con  el  cambio  y  —  ux  podemos  poner  y'  =  2-  —  (- )2  —  2 u  —  u2 .  Como  y'  —  u'x  +  u. 
sustituyendo  tenemos  u'x  +  u  —  2u  —  u2 ,  es  decir,  xu'  —  u  —  u2 . 

Si  u  ^  u2,  podemos  poner  Para  integrar,  descomponemos  — j~2  =  ^  + 

lo  que  se  satisface  para  A  =  B  —  1.  Entonces,  integrando,  logw  —  log(l  —  u)  =  log  fj.  es 

decir,  y  sustituyendo  u  =  ^  tenemos  x-y/x  ~  ü’  de  donde  Cy  —  x(x  —  y).  De 

aquí  es  fácil  despejar  explícitamente  y  si  así  se  desea. 

Por  otra  parte,  a  partir  de  Uo  =  0  y  Uq  =  1  (para  las  cuales  u  —  u 2),  se  tienen  las 
soluciones  singulares  y  —  0  e  y  —  x. 


APARTADO  3'. 

Reducibles  a  homogéneas. 

Consideremos  la  ecuación 


t  _  f  ( a^x  +  &i y  + 

\  ax  +  by  +  c  ) 

Para  resolverla,  hay  que  distinguir  dos  casos: 

3'.1.  Supongamos  en  primer  lugar  que  las  rectas  ax  +  by  +  c  =  0  y  a\x  +  b\y  +  c\  —  0 
se  cortan  en  el  punto  (xo,  yo)-  Así,  tendremos  que  ax  +  by  +  c  =  a(x  —  xo)  +  b(y  —  yo)  y 
a\x  +  biy  +  c\  —  ai (x  —  xo)  +  b± (y  —  yo) .  Hagamos  ahora  el  cambio  de  variable  y  de  función 
X  =  x  —  Xo,  Y  —  y  —  yo,  con  lo  cual 


y,  _  ,  _  r  ( ai(x  -  a^o)  +  bi(y  -  y0)\  _  / atX  +  hY\  =  í ai+btY\ 

J  V  a(x  -  ®o)  +  b(y  -yo)  )  \  aX  +  bY  )  y  a  +  bx  J  1 

es  decir,  hemos  reducido  la  ecuación  a  una  homogénea. 

3'. 2.  En  segundo  lugar,  supongamos  que  ax  +  by+c  —  0  y  aix  +  biy  +  ci  =  0  son  rectas 

paralelas,  con  lo  cual  podrá  ponerse  (ai,  b\)  —  K(a,  b)  para  algún  K  e  R.  Efectuemos  ahora 

/ _ 

el  cambio  de  función  z  —  ax  +  by.  Derivando,  z  =  a  +  by  ,  o  sea,  y  =  Si  sustituimos 

en  la  E.  D.  original  obtenemos 


dz 

dx 


a  +  bf 


Kz  +  ci 


que  es  de  variables  separadas. 


z  +  c 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


RECETA  3'.  Reducibles  a  homogéneas. 

Son  de  la  forma 

,  _  ,  í a\x  +  biy  +  ci\ 

\  ax  +  by  +  c  J 

3'.1.  Si  las  rectas  ax  +  by  +  c  —  0  y  a\x  +  b\y  +  c\  —  0  se  cortan 
en  ( xo,yo ),  se  hace  el  cambio  de  variable  y  de  función  X  —  x  —  xq, 
Y  —  y  —  yo-  La  ecuación  se  reduce  a  una  homogénea. 

3'. 2.  Si  ax  +  by  +  c  —  0  y  aix  +  biy  +  ci  —  0  son  rectas  paralelas,  se 
hace  el  cambio  de  función  z  —  ax  +  by.  La  nueva  ecuación  que  aparece 
es  de  variables  separadas. 


Ejemplo  3'.1.  Resolver 


— 2x  +  4y  —  Q 
x  +  y  —  3 


Las  rectas  —2x  +  Ay  —  6  —  0  y  x  +  y  —  3  =  0  se  cortan  en  el  punto  (x,  y)  =  (1,  2),  con 

lo  que  efectuamos  el  cambio  X  —  x  —  1,  Y  =  y  —  2,  Y'  =  y' .  Sustituyendo,  obtenemos  la 

ecuación  homogénea 

,  _  -2X  +  4 Y  _  -2  +  4 Y/X 
~  X  +  Y  ~  1  +  Y/X 

Para  resolverla,  hacemos  un  nuevo  cambio  u  —  Y/X ,  de  donde  Y  =  uX,  Y'  =  u'X  +  u. 
Tras  sustituir,  tenemos  u'X  +  u  —  ~i^u  que,  a  la  postre,  podemos  poner  como  —X  ^  = 

u2 -3u- \-2 
u-\- 1 

Tenemos  ahora  que  distinguir  cuándo,  y  cuándo  no,  se  anula  la  expresión  u2  —  3u  +  2. 
Resolviendo  u2  —  3u  +  2  =  0,  esto  ocurre  para  u  —  1  y  u  —  2. 

Analicemos  en  primer  lugar  el  caso  u2  —  3u  +  2  ^  0.  Así,  podemos  escribir 


dX 

~Y 


—  (u  +  1)  du 
u2  —  3u  +  2 


A 

u-1 


du  + 


B 

u  —  2 


du  = 


2 


u-1 


du  — 


3 

u  —  2 


du 


de  donde,  integrando,  21og(,u  —  1)  —  31og(it  —  2)  =  log(JiA)  y,  consiguientemente, 
=  KX.  Sustituyendo  ahora  u  =  Y/X  llegamos  fácilmente  a  (Y— X)2  =  K(Y— 2A)3; 
y  volviendo  a  las  variables  originales  x  e  y  obtenemos  las  soluciones  (y—x—1)2  —  K(y— 2x)3 
de  la  E.  D.  de  partida. 

Finalmente,  con  uo  =  1  y  uq  —  2  tenemos,  respectivamente,  las  soluciones  Y  =  X  e 
Y  —  2X  que,  sustituyendo  X  e  Y  por  su  valor,  se  traducen  eny  —  x  +  lyy  —  2x. 


E.  D.  explícitas  de  primer  orden 


11 


Ejemplo  3'. 2.  Resolver 

/  x-y-1 

y  = - o- 

x  —  y  —  2 


Efectuamos  el  cambio  de  función  z  =  x  —  y,  de  donde  z'  =  1  —  y' .  Sustituyendo, 
tenemos  —z'  +  1  =  zE^,  Y  consiguientemente  =  —  1  =  z^,  que  es  la  ecuación 

(z  —  2)  dz  =  —dx,  cuyas  variables  están  separadas.  Integrando,  \(z  —  2)2  =  —  x  +  K,  y 
finalmente,  sustituyendo  de  nuevo  z  —  x  —  y  y  denotando  C  —  2 K,  obtenemos  que  las 
soluciones  de  la  E.  D.  original  son  (x  —  y  —  2)2  +  2x  =  C. 


APARTADO  3". 


Homogéneas  implícitas. 

Sea  la  ecuación 


Para  resolverla,  consideremos  la  curva  F(a,f3 )  =  0  y  supongamos  que  hemos  logrado 
encontrar  una  representación  paramétrica  de  la  curva  dada  por  a  =  <p(t),  ¡3  —  p(í).  Es 
decir,  que  se  verifica  F(<p(t),tp(t))  =  0.  Hagamos  ahora  el  cambio  de  función  y  por  t 
mediante  |  =  <p(t),  teniendo  en  cuenta  que  y'  =  ^(t). 

Si  derivamos  y  =  x(p(t)  respecto  de  x  tenemos  y'  —  p(t)  +  x(p'{t)jr¿,  es  decir, 
V>(í)  =  <p(t)  +  o 

^(í)  -  <p(t)  =  x<p'(t)^ 

que,  en  principio,  es  una  ecuación  en  variables  separadas. 

•  Si  ip(t)  <p(t),  podemos  poner  cuya  solución  será  x  =  con 

(f)(t)  —  f  De  aquí  hue  D-  D.  de  partida  tiene  las  soluciones 


x  =  Cem 
y  =  C(p{t)e m  ’ 


•  Si  existe  t0  tal  que  ip(to)  —  <p(t o),  formalmente,  podemos  pensar  0  =  x(p'(t)^,  luego 
<p'(t)  —  0  y  por  tanto  ip{t)  —  cte.  =  <p(to),  lo  que  nos  llevaría  a  la  solución  y  —  x<p(to).  Esta 
recta  es,  efectivamente,  una  solución  de  la  E.  D.,  como  podemos  comprobar  directamente: 
F(l,y')  =  F(<P(to),<p(t0))  =  F(ip(t0),ip(t0))  =  0. 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


RECETA  3".  Homogéneas  implícitas. 

Son  de  la  forma 


Consideramos  la  curva  F(a,  ¡3)  =  0.  Si  encontramos  una  representación 
paramétrica  a  —  <p(t),  (3  —  ^(t),  F(<p(t),ip(t))  =  0,  se  hace  el  cambio  de 
función  y  por  t  mediante  |  y'  =  ^(t).  Así,  derivando  y  —  x<p(t ) 

respecto  de  x,  aparece  una  ecuación  en  variables  separadas. 


Ejemplo  3".  Resolver 

x2(y')2-(y2+x2)=  0. 

Si  ponemos  la  ecuación  en  la  forma  (y')2  —  (-)2  =  1,  podemos  recordar  que  el  coseno 
y  el  seno  hiperbólicos  satisfacen  la  relación  (chí)2  —  (shí)2  =  1,  lo  que  se  adecúa  a  nuestras 
necesidades.  Así,  tomemos  ahora  -  =  shí  y  y'  =  chf.  Si  derivamos  y  —  xsht.  tenemos 
y ’  —  shí  +  xcht  4¿,  o  sea,  chí  =  sh t  +  x ch t  ^ .  Despejando,  ^  =  ch t  dt  (notar  que, 
en  esta  ecuación,  el  denominador  no  se  anula  nunca  ya  que  chí  >  sh  t.  luego  no  hay  que 
preocuparse  de  analizar  por  separado  las  raíces  de  ch  t  —  sh  t  —  0)  e,  integrando,  x  —  Ce 
con 

m  =  /  d nvW  dt  =  ¡  e-i^- dt  =  5  /(1  +  e2<)  dt  =  h  +  h2‘- 

Entonces,  la  E.  D.  original  tiene  como  soluciones  las  curvas 

Jx  =  C'eí/2+e2‘/4 
\  y  =  C(shí)eí/2+e2í/4 


APARTADO  3"'. 

Sea  la  ecuación  y ’  =  f(x,y ) 

con  /  tal  que,  para  algún  ct  fijo,  verifica 

f(Xx,Xay)  =  X  a_1/(x,y). 

Nótese  en  primer  lugar  que,  cuando  a  —  0,  sin  más  que  tomar  X  —  x  tenemos  y'  =  f(x,  y)  = 
x~x  f(l,y),  que  es  una  ecuación  en  variables  separadas;  y,  cuando  a  —  1,  se  puede  poner 

y'  =  f{x,y )  =  /  (^D  =  /  (i»  |)  • 

es  decir,  nos  encontramos  ante  una  E.  D.  homogénea.  En  otro  caso,  veamos  cómo  el  cambio 
de  función  y  —  za  transforma  la  ecuación  en  una  homogénea: 


E.  D.  explícitas  de  primer  orden 
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Derivando,  y'  =  aza  lz'  y,  sustituyendo  en  la  E.  D.  original,  aza  lz'  —  f(x,za),  es 


decir, 


¿  =  -(-)Q_1/(x,^)  =  -f  í-x,  (-) >)  =  -/  (-,l) 
ot\zJ  a  \z  \zJ  J  a  \z  / 


que,  efectivamente,  es  una  E.  D.  homogénea. 

Lógicamente,  como,  al  hacer  en  la  homogénea  el  cambio  z  =  ux,  ésta  se  transforma  en 
una  de  variables  separadas,  si  hubiéramos  efectuado  desde  el  principio  el  cambio  y  =  ( ux)a , 
nuestra  E.  D.  se  hubiera  convertido  directamente  en  una  de  variables  separadas. 


Por  último  comentar  que,  extrayendo  X  —  x,  este  tipo  de  ecuaciones  puede  ponerse 
como 

y '  =  =  f  JO  =  xa~xf  (J  JO  =  xa~l]l  (JO  • 

Pero,  al  menos  a  simple  vista,  no  parece  más  sencillo  describir  las  ecuaciones  que  estamos 
tratando  como  “las  que  tienen  la  forma  y'  =  xa~xh{yx~a)n  que  tal  como  lo  hemos  hecho 
en  este  apartado. 


RECETA  3"'.  Si  la  ecuación  y'  —  f(x,y) 
es  tal  que,  para  algún  a  ^  0  fijo,  /  satisface 

f(Xx,Xay)  =  X  a~xf(x,y), 

entonces  el  cambio  de  función  y  —  za  transforma  la  ecuación  en  una 
homogénea.  (Si  ct  =  1,  la  E.  D.  ya  es  homogénea;  y  si  /  cumple  la 
relación  anterior  con  a  —  0,  la  E.  D.  es  de  variables  separadas.) 


Ejemplo  3"'.  Resolver 


/ 

y 


1  y  _  0  y/x 

2  x  y2 


Dada  /(x,y)  =  \x  xy  —  3xl/2y  2,  para  intentar  encontrar  a  tanteamos  /( Xx,  X ay)  = 
^(Xx)~1(Xay)  —  3(Xx)1/2(X°íy)-‘2  —  \Xa~1x~1y  —  ?>X1/2~2oix1/2y~2 ,  y  observamos  que  esto 
es  igual  a  A a~1f(x,y)  sin  más  que  tomar  a  — 

Entonces,  si  sustituimos  y  =  z1!2  en  la  E.  D.,  tenemos  \z~x!2z'  =  \x~xz1^2  — 
2>xx^2z~x,  es  decir,  z'  —  |  —  G(^)^1/2,  que  es  homogénea.  Para  resolverla,  tomamos  ahora 
z  =  ux ,  con  lo  cual,  sustituyendo,  u'x+u  —  u  —  6u~x^2:  o  sea,  u 1/2  du  =  —  &x~x  dx,  ecuación 
en  variables  separadas.  Integrándola,  | v?!2  =  — 61og(^),  luego  x  =  Cexp(— -u3/2/9). 
Deshaciendo  los  cambios  z  =  ux  y  y  —  ux  encontramos  que  las  soluciones  de  la  E.  D. 
de  partida  son,  expresadas  como  curvas  en  paramétricas, 

íx  =  Ce~u3/2/9 

{  y  =  G i/V/V«s/2/i8  ' 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


APARTADO  4. 

Ecuaciones  exactas. 

Llamamos  exacta  a  una  ecuación  diferencial 

P(x,  y)  dx  +  Q(x,  y)  dy  =  0, 

es  decir,  y'  =  ^  ,  que  cumple  Py  =  Qx  (con  la  notación  Py  =  ,  Qx  =  §§). 

Antes  de  explicar  como  resolverlas,  comentemos  brevemente  algo  sobre  “expresiones 
diferenciales”  (rigurosamente  hablando,  estamos  tratando  con  1  -formas  diferenciales  w  = 
P  dx  +  Q  dy,  aunque  no  entraremos  en  ello): 

Supongamos  de  antemano  en  todo  lo  que  sigue  que  P  y  Q  son  de  clase  C1  (continuas 
con  derivadas  parciales  continuas)  en  su  dominio  de  definición  (un  abierto  de  M2).  Una 
expresión  diferencial  P(x,  y)  dx  +  Q{x,  y)  dy  se  dice  que  es  una  diferencial  cerrada  en  una 
región  R  del  plano  xy  si  se  verifica  Py(x,y )  =  Qx(x,y )  para  todo  ( x,y )  G  R.  Y  se  dice 
exacta  en  R  cuando  existe  alguna  función  F(x,y )  tal  que  ^  =  P  y  ^  =  Q  para  todo 
(x,  y)  G  R;  en  otras  palabras,  si  la  diferencial  de  F  es  dF  —  P  dx  +  Q  dy  (F,  que  es 
única  salvo  constantes,  se  denomina  función  potencial).  El  teorema  de  Schwartz  sobre 
igualdad  de  derivadas  cruzadas  nos  asegura  que  cualquier  expresión  diferencial  exacta 
es  cerrada.  Lo  contrario  no  es  cierto  en  general,  aunque  sí  en  una  clase  muy  amplia 
de  dominios  de  IR2:  los  simplemente  conexos  que,  intuitivamente,  son  los  que  no  tienen 
agujeros.  Demostrar  este  hecho  uo  es  excesivamente  sencillo,  pero  tampoco  es  necesario 
para  lo  que  aquí  pretendemos.  En  realidad,  el  lema  de  Poincaré  (que  normalmente  se 
prueba  en  cualquier  curso  de  cálculo  integral  en  varias  variables)  asegura  que  una  expresión 
cerrada  es  exacta  siempre  que  el  dominio  sea  estrellado ,  lo  que  significa  que  exista  un 
punto  del  dominio  que  se  pueda  unir  a  todos  los  demás  mediante  un  segmento  sin  salimos 
del  dominio;  en  particular,  los  conjuntos  convexos  son  estrellados.  Además,  esto  asegura 
que,  dada  cualquier  expresión  cerrada,  es  exacta  localmente,  es  decir,  alrededor  de  cada 
punto  podemos  restringir  el  dominio  de  tal  forma  que  la  expresión  sea  exacta  en  ese  nuevo 
dominio.  Por  lo  tanto,  en  lo  que  a  nosotros  concierne,  podemos  identificar  los  conceptos 
de  exacto  y  cerrado,  ya  que  no  nos  estamos  preocupando  de  dónde  están  definidas  las 

E.  D.  que  tratamos  de  resolver  ni  en  qué  intervalo  existen  las  soluciones.  En  realidad,  en 
ecuaciones  diferenciales  suele  hablarse  siempre  de  exacto  aún  refiriéndose  a  que  se  satisface 
la  igualdad  Py  =  Qx. 

Una  E.  D.  exacta  es  una  expresión  exacta  igualada  a  cero.  Veamos  como  resolverlas: 
si  tenemos  P  dx  +  Q  dy  =  0  exacta,  como  existe  F  tal  que  dF  =  P  dx  +  Q  dy,  entonces  la 
ecuación  podemos  ponerla  en  la  forma  dF  =  0  y,  por  tanto,  su  solución  será  F(x,y)  =  C 
(siendo  C  constante  arbitraria).  Así  pues,  basta  con  que  encontremos  la  función  potencial 

F .  El  procedimiento  para  hallarla  que,  tal  como  veremos,  funciona  gracias  a  que  Py  —  Qx, 
es  como  sigue: 

Buscamos  F  tal  que  jP  —  P:  así,  es  posible  encontrar  F  integrando  P[x,  y)  respecto 
a  x  mientras  se  mantiene  y  constante,  es  decir,  F(x,y)  —  f  P(x,y)dx  +  <p(y),  donde  la 
función  arbitraria  ip(y)  es  la  “constante”  de  integración.  Derivando  respecto  de  y  obtenemos 


E.  D.  explícitas  de  primer  orden 
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=  Ey  f  P(X’V)  d'x  +  ‘P'i.v)-  P°r  °U’a  parte,  si  utilizamos  ^  =  Q,  de  aquí  resulta  que 
¡p'( y )  =  Q(x,y)  ~  -§¡j  J  P{x,y)  dx ;  ésta  es  realmente  una  expresión  independiente  de  x  ya 
que 


d_ 

dx 


d_ 

dy 


dQ  d  /  d 
dx  dy  V  dx 


dQ  dP 
dx  dy 


Una  vez  conocida  < p'(y ),  integrando  obtenemos  <p(y)  y,  sustituyendo  su  valor,  llegamos  a 
la  función  potencial  F(x,y).  Así,  quedan  halladas  completamente  las  soluciones  buscadas 
F(x,y )  =  C,  expresadas  en  forma  implícita. 

(Lógicamente,  para  encontrar  F  podría  haberse  seguido  el  proceso  anterior  cambiando 
el  orden  en  el  que  se  usa  P  y  Q,  partiendo  de  =  Q  y  posteriormente  utilizando  ^  =  P. 
Asimismo,  integrando  estas  dos  expresiones  e  igualándolas,  muchas  veces  basta  una  simple 
inspección  para  determinar  F .) 


RECETA  4.  Ecuaciones  exactas. 

Son  las  de  la  forma 

P(x,  y)  dx  +  Q(x,  y)  dy  —  0, 

es  decir,  y'  —  ¿y  —  —  Qpé’y]  ’  due  cumplen  Py  —  Qx.  Se  busca  una 
función  F(x,  y)  tal  que  dF  =  u>  =  P  dx  +  Q  dy,  y  la  solución  de  la  E.  D. 
es  F(x,y )  =  C  (siendo  C  constante). 


Ejemplo  4.  Resolver 


3y  +  ex  +  (3x  +  eos  y)y'  —  0. 


Si  ponemos  la  ecuación  en  la  forma  P  dx  F  Q  dy  =  0  con  P(x,y)  =  3y  +  ex  y 
Q(x,y)  —  3x  +  eos  y,  es  claro  que  Py  —  Qx  =  3,  luego  la  E.  D.  es  exacta.  Calculemos 
la  función  potencial  F  (que  nos  dará  directamente  las  soluciones  F(x,y )  =  C).  Como 
Fx  =  3y  +  ex,  integrando  respecto  de  x,  F(x,  y)  —  3 yx  +  ex  +  <p(y).  Derivando  respecto  de 
y  e  igualando  a  Q  queda  3x  +  <p'(y)  —  3x  +  eos  y,  es  decir,  ¡p'(y)  —  eos  y,  de  donde  basta 
tomar  <p(y)  =  sen  y,  y  por  tanto  F(x,y )  =  3 yx  +  ex  +  sen  y.  Así,  la  solución  de  la  E.  D. 
viene  dada,  implícitamente,  por  3yx  +  ex  +  sen  y  —  C.  (Nótese  que  al  integrar  <p'(y)  —  eos  y 
no  hace  falta  poner  la  constante  de  integración  íp(y)  —  sen  y  +  C\  ya  que,  en  ese  caso,  una 
de  las  dos  constantes  de  la  solución  3 yx  +  ex  +  sen  y  +  C\  —  C  sería  claramente  superflua.) 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


APARTADO  4'. 

Reducibles  a  exactas:  Factores  integrantes. 

Si  tenemos  una  ecuación  P(x,  y)  dx  +  Q(x,  y)  dy  —  0  que  no  es  exacta,  una  idea  para 
intentar  resolverla  sería  tratar  de  encontrar  alguna  función  y(x,y)  no  idénticamente  nula 
tal  que 

y(x,  y)P(x,  y)  dx  +  y(x,  y)Q(x,  y)  dy  =  0 

sea  exacta.  Como  esta  ecuación  es  equivalente  a  la  de  partida,  sus  soluciones  y  las  de 
P(x,  y)  dx  +  Q(x,  y)  dy  —  0  serán  las  mismas. 

Desgraciadamente,  no  hay  ningún  procedimiento  general  que  permita  encontrar 
factores  integrantes.  Sin  embargo,  sí  que  es  posible  hacerlo,  y  de  manera  sencilla,  en  dos 
casos: 

4'.1.  Existencia  de  factor  integrante  de  la  forma  y(x).  Queremos  que  / jl{x)P{x ,  y)  dx  + 
y(x)Q(x,y )  dy  —  0  sea  exacta,  esto  es, 

d  d 

—  (y(x)P{x,y))  =  —  (y(x)Q(x,y)). 

Derivando,  y(x)Py(x,y)  =  y'(x)Q(x,  y)  +  y(x)Qx(x,  y),  osea,  y(x)(Py(x,  y)  -  Qx(x,  y))  = 
y'(x)Q(x,y).  Para  que  esto  tenga  sentido, 

h'O)  =  Py(x,y)  -  Qxjx,  y) 
y(x)  Q{x,y ) 

tiene  que  resultar  ser  una  función  que  dependa  exclusivamente  de  x,  que  denotamos  h(x). 
Cuando  éste  es  el  caso,  es  claro  que  la  función  y  que  satisface  la  relación  anterior  es 

y(x)  —  exp  (/  h{x)dx 

con  lo  cual  hemos  encontrado  el  factor  integrante  buscado. 

4'.  2.  Existencia  de  factor  integrante  de  la  forma  y  (y).  Repitiendo  el  proceso  anterior 
buscamos  -§¿(y(y)P(x,y))  =  £(y(y)Q(x,y)),  es  decir,  y'(y)P  +  y(y)Py  =  y(y)Qx ,  Y  por 

tanto  =  Q* pPy ,  que  tiene  que  ser  función  sólo  de  y,  que  denotamos  h(y).  En  estas 
condiciones,  el  factor  integrante  es  y(y)  =  exp (f  h(y)  dy). 

4'. 3.  Aparte  de  los  casos  anteriormente  tratados,  para  algunos  tipos  de  problemas  se 
puede  intentar  encontrar  factores  integrantes  imponiendo  a  y(x,y)  condiciones  restrictivas 
de  muy  diverso  tipo.  Por  ejemplo,  exigiendo  que  sea  de  la  forma  y(x,y)  —  xay13  con  a 
y  ¡3  constantes  a  determinar,  que  sea  y(x  +  y),  o  y(xy),  etc.  Para  estudiar  estos  casos, 
lo  que  hay  que  hacer  siempre  es  igualar  J^(yP)  =  yQ )  e  intentar  resolver  la  nueva 
ecuación  que  aparece,  teniendo  en  cuenta,  sobre  todo,  si  tiene  sentido.  Por  ser  generalmente 
procedimientos  bastante  particulares,  no  comentaremos  aquí  nada  más  sobre  ellos. 


E.  D.  explícitas  de  primer  orden 


17 


RECETA  4'.  Reducibles  a  exactas:  Factores  integrantes. 

Si  P(x,y)dx  +  Q(x,y)dy  =  0  no  es  exacta,  podemos  intentar 
encontrar  p(x,y)  tal  que 

p(x,  y)P{x,  y)  dx  +  p(x,  y)Q{x,  y)dy  =  0 

sea  exacta. 

4,.l.  Existencia  de  factor  integrante  de  la  forma  p(x).  Ocurre 
cuando  Py  =  h(x),  tomándose  p{x)  —  exp (J  h(x)  dx). 

4'. 2.  Existencia  de  factor  integrante  de  la  forma  p(y).  Ocurre 
cuando  ^pPg  =  h(y),  tomándose  p{y)  =  exp (J  h(y )  dy). 

4'. 3.  Otras  expresiones  restrictivas  para  p(x,y). 


Ejemplo  4'.  Resolver 

( 2x 2  +  y)  dx  +  ( x2y  —  x)  dy  —  0. 

En  este  caso,  P(x,y)  —  2x2  +  y  y  Q(x,y)  —  x2y  —  x.  Esta  ecuación  no  es  exacta  ya 
que  Py  =  1  y  Qx  —  2 xy  —  1.  Para  intentar  encontrar  un  factor  integrante  se  calcula 

Py  ~  Qx  _  1  -  (2 xy  -  1)  _  2(1  -  xy)  _  -2 
Q  x2y  —  x  —x(l—xy)  x 

Ya  que  se  obtiene  una  expresión  que  depende  sólo  de  x,  podemos  asegurar  que  existe 
un  factor  integrante  dado  por  la  fórmula  p(x)  =  exp (J  ^  dx)  —  x~2 .  Entonces,  si 
multiplicamos  la  E.  D.  por  p(x)  =  x~2  se  obtiene  la  ecuación  exacta  (2  +  yx~2)  dx  +  (y  — 
x~x)  dy  —  0.  Por  el  método  usual,  encontramos  que  la  función  F  tal  que  Fx  —  2  +  yx~2  y 
Fy  =  y  —  x~l  es  F(x,y)  =  2x  —  yx~x  +  \y2 ■  Por  tanto,  la  solución  de  la  E.  D.  exacta,  y 
también  de  la  de  partida,  resulta  ser  2x  —  yx~l  +  \y2  —  C. 

APARTADO  5. 

Ecuaciones  lineales  de  primer  orden. 

Dada  la  ecuación 

y'  +  a(x)y  =  b{x), 

vamos  a  explicar  cómo  resolverla  por  tres  métodos  distintos: 

(i)  Encontrar  un  factor  integrante  de  la  forma  p(x).  Para  ello,  si  la  ponemos  en  la 
forma  ( a(x)y  —  b(x))  dx  +  dy  =  0  y  denotamos  P(x,y)  —  a(x)y  —  b(x)  y  Q(x,y)  —  1,  se 
tiene  v qx  =  a(x).  Por  tanto,  según  liemos  visto  anteriormente,  la  E.  D.  tiene  el  factor 
integrante  p(x)  —  exp (f  a(x)dx).  Así,  multiplicando  por  p(x),  la  ecuación 

a(x)  dx  \  ( a(x)y  —  b(x))  dx  +  exp 


exp 


a(x)  dx  \  dy  —  0 
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tiene  que  ser  exacta.  Ahora,  bastará  encontrar  la  función  potencial  F  con  lo  que  la  ecuación 
anterior  podrá  ponerse  dF  =  0  y  su  solución  será  F(x,  y)  =  C.  Busquemos  F: 


Como  Fy  =  exp (f  a(x)dx),  tendremos  F  =  yexp(f  a(x)  dx)  +  <p(x).  Por  otra  parte, 
derivando  esta  F  respecto  de  a;  y  usando  que  Fx  —  exp (f  a(x)  dx)(a(x)y  —  b(x ))  llegamos 
a 


Fx  =  y  exp  (  /  a(x)  dx  )  a(x)  +  <p'(x)  —  exp  (  /  a(x)  dx  )  ( a(x)y  —  b(x)), 


de  donde  —b(x)exp(f  a(x)  dx)  —  <p'(x).  Integrando,  <p(x)  —  —fb(x)exp(f  a(x)  dx)  dx, 
luego  F  =  yexp(f  a(x)  dx)  —  f  b(x)  exp (f  a(x)  dx)  dx  y  la  solución  de  la  ecuación  exacta 
(y  de  la  lineal  de  partida)  es,  expresada  en  forma  implícita, 


y  exp  (  /  a(x)  dx)  —  b(x)  exp  (  /  a(x)  dx)  dx  —  C. 


Sin  más  que  despejar  y ,  tenemos  que  la  solución  de  la  ecuación  lineal  resulta  ser 


y  =  exp 


a(x)  dx 


j b(x) exp 


a(x)  dx 


dx  +  C 


(ü)  Un  segundo  método  de  resolución  se  basa  en  resolver  previamente  la  ecuación 
lineal  homogénea  asociada  y'  +  a(x)y  —  0.  Esta  ecuación  es  de  variables  separadas,  pues 
puede  ponerse  ^  =  —a(x)  dx;  su  solución  es  y  —  (7exp(—  f  a(x)  dx). 

Apliquemos  ahora  el  método  de  variación  de  las  constantes ,  esto  es,  consideremos 


y  —  C(x)  exp 


a(x)  dx 


y  vamos  a  ver  cómo  debe  ser  C(x)  para  que  se  verifique  y'  +  a(x)y  —  b(x).  Derivando, 
y'  =  C'{x)  exp(—  f  a(x)  dx)  —  C(x)a(x)  exp(—  f  a(x)  dx)  y,  sustituyendo  en  la  ecuación 
lineal, 


C\x)  exp  (  —  /  a(x)  dx  )  —  C(x)a(x)  exp 


a(x)  dx 


+  a(x)C(x)  exp 


a(x)  dx 


b(x). 


Dos  de  los  sumandos  anteriores  se  cancelan,  de  donde  C'{x)  =  b(x)exp(f  a(x)  dx)  e, 
integrando, 

C(x)  —  í  b(x)  exp  f  í  a(x)  +  C. 


Así,  hemos  llegado  a  la  misma  expresión  para  las  soluciones  que  la  que  encontramos  por 
el  método  anterior. 
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(iii)  El  tercer  procedimiento  de  resolución  parte  de  suponer  que  liemos  encontrado, 
de  alguna  forma,  una  solución  particular  yp(x)  de  la  E.  D.  lineal.  Entonces,  la  solución 
general  de  la  lineal  es  yp  más  la  solución  general  de  la  lineal  homogénea  asociada,  es  decir, 


y  =  yv  +  C  exp 


a(x)  dx 


es  solución  para  todo  C  e  R.  La  justificación  de  este  hecho  es  sencilla.  En  efecto,  basta 
comprobar  que,  si  yp  es  solución  de  y'  +  a(x)y  —  b( x)  y  y  lo  es  de  y'  +  a(x)y  —  0,  entonces 
y  +  yp  es  solución  de  y'  +  a(x)y  —  b(x),  lo  cual  es  claramente  cierto: 

(■ y  +  VpY  +  a(x  )(y  +  yp)  =  (; y '  +  a(x)y)  +  (yp  +  a(x)yp)  =  0  +  b(x)  =  b(x). 

(iv)  El  último  método  de  resolución  de  ecuaciones  lineales  que  describimos  consiste 
en  efectuar  una  descomposición  y(x)  —  u(x)v(x)  adecuada.  Tomando  y  de  esa  forma,  si 
derivamos,  y'  =  u'v  +  uv' ,  con  lo  cual,  al  sustituir  en  la  ecuación,  u' v  +  uv'  +  a(x)uv  =  b(x). 
Sacando  u  factor  común,  podemos  escribir  la  expresión  anterior  como  u' v  +  (v'  +  a(x)v)u  = 
b(x).  Vamos  ahora  a  elegir  v  de  tal  forma  que  se  anule  el  coeficiente  de  u,  es  decir, 
que  satisfaga  v’  +  a(x)v  =  0.  Esta  es  una  E.  D.  en  variables  separadas;  resolviéndola, 
^  =  —a(x),  con  lo  cual  basta  tomar  logu  =  —  f  a(x)  dx,  es  decir, 


v(x)  —  exp 


a(x)  dx 


Con  v(x)  esa  función,  la  ecuación  queda  ahora  u'v  =  b(x),  de  donde  u'  —  b(x)v  l,  es  decir, 
u'(x)  =  b(x)exp(f  a(x)dx).  Integrando, 


u(x) 


b(x)  exp 


a(x)  dx 


dx  +  C. 


Sin  más  que  recomponer  y  =  u(x)v(x),  con  este  procedimiento  de  nuevo  encontramos  la 
misma  expresión  para  las  soluciones  de  la  E.  D.  lineal  de  partida. 

Para  concluir,  comentemos  una  vez  más  que  no  hace  falta  recordar  la  fórmula  que 
hemos  obtenido  para  las  soluciones  de  la  ecuación  lineal,  sino  que  basta  seguir  en  cada 
problema  alguno  de  los  procesos  descritos.  Generalmente,  en  opinión  del  que  escribe,  los 
que  suelen  conducir  a  la  solución  por  un  procedimiento  más  corto  suelen  ser  el  segundo  y 
el  cuarto.  El  tercero  tiene,  sobre  todo,  gran  importancia  teórica.  Además,  merece  la  pena 
destacar  que  el  segundo  y  el  tercero  tienen  su  paralelismo  a  la  hora  de  resolver  ecuaciones 
lineales  de  orden  superior  (como  veremos  más  adelante) ,  mientras  que  los  otros  dos  sólo  se 
aplican  a  las  de  primer  orden. 

Ejemplo  5.  Resolver 


por  los  cuatro  métodos  descritos. 


2  xy'  —  3  y  =  4x2 
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RECETA  5.  Ecuaciones  lineales  de  primer  orden. 

Son  de  la  forma 

y'  +  a(x)y  =  b(x). 

Hay  tres  métodos  de  resolución:  (i)  Encontrar  un  factor  integrante 
de  la  forma  y(x).  (■ ii )  Resolver  la  ecuación  lineal  homogénea  asociada 
y'  +  a(x)y  —  0  (que  es  de  variables  separadas),  cuya  solución  es 
y  =  Cexp(—  f  a(x)  dx ),  y  usar  el  método  de  variación  de  las  constantes 
(esto  es,  cambiar  C  por  C(x)  en  la  expresión  anterior  y  sustituir 
en  la  ecuación  lineal),  (iii)  Encontrar  de  alguna  forma  una  solución 
particular  yp(x),  con  lo  cual  la  solución  general  de  la  lineal  es  yp 
más  la  solución  general  de  la  homogénea  asociada,  (iv)  Descomponer 
y(x)  =  u(x)v(x),  sustituir  en  la  lineal,  e  igualar  a  0  el  coeficiente  de 
u ,  resolviendo  la  ecuación  que  aparece  (vr  +  a(x)v  =  0,  que  es  de 
variables  separadas);  tras  esto,  queda  una  ecuación  en  u(x)  de  variables 
separadas. 

De  cualquier  modo  se  obtiene  que  la  solución  general  de  la  E.  D. 
lineal  es 


y  =  exp 


a(x)  dx 


b(x)  exp  (  /  a(x )  dx)  dx  +  C 


(i)  Tenemos  la  ecuación  lineal  y'  —  y  =  2x,  es  decir,  dy  —  (j-y  +  2 x)  dx,  que 
es  de  la  forma  P  dx  +  Q  dy  =  0  con  P(x,y)  =  p^y  —  2x  y  Q(x,y)  =  1.  Como 

Py  q®x  =  3^2^  0  =  existe  el  factor  integrante  y(x)  =  exp (f  dx)  —  x~3//2.  Así,  la 
ecuación  (=p x~r^2y  —  2x~ l/2)  dx  +  x~ 3/2  dy  =  0  es  exacta.  La  función  potencial  F  debe 
cumplir  Fy  —  a:-3/2,  luego  F  —  x~3/,2y  +  ip(x).  Por  otra  parte,  Fx  =  =^-x~5^2y  —  2x~1^2  = 
=Y-x~5^2y  +  tp'(x),  de  donde  ip(x)  =  —2  f  x~x!2  dx  =  —4a:1/2.  Por  tanto,  F(x,y)  — 
x~3/2y  —  4a:1//2,  y  la  solución  de  la  E.  D.  es  x~3!2y  —  4a:1/2  —  C,  o  sea  y  —  Cx 3/2  +  4a:2, 
del. 


(ii)  La  lineal  homogénea  asociada  es  y'  —  ^ y  —  0.  Podemos  ponerla  como  ^ 

de  variables  separadas,  cuya  solución  es  log y  =  |loga:  +  K ,  es  decir,  y  —  Cx 3/2. 
Empleemos  ahora  el  método  de  variación  de  las  constantes,  para  lo  cual  tomamos  y  = 
C(x)x3/2 .  Derivando,  y'  =  C' (x)x3¡2  +  3C(x)xx /2  y,  sustituyendo  en  E.  D.  de  partida, 
2 x(C' (x)x3t2  +  3C(x)x 3/2)  —  3C(x)x3/2  —  4a:2,  esto  es,  C'(x)  =  2a:-1/2.  Integrando, 
C(x)  =  4a:1/2  +  Ki  luego  la  solución  de  la  lineal  es  y  =  (4a:1/2  +  K i)a:3/2.  Si  empleamos 
de  nuevo  C  para  denotar  la  constante,  y  =  Cx 3/2  +  4a:2,  la  misma  expresión  que  hemos 
encontrado  en  (i). 


(iii)  Primero,  tratemos  de  hallar  una  solución  particular  de  la  ecuación  lineal.  Lo  más 
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sencillo  es  intentar  probar  si  existe  alguna  solución  polinómica.  Esto  sólo  sería  posible 
con  un  polinomio  de  segundo  grado,  pues  en  otro  caso  no  podrían  cancelarse  nunca  todos 
los  sumandos.  Por  tanto,  vamos  a  tantear  con  polinomios  de  la  forma  y  =  ax2  +  bx  +  c. 
Derivando,  y'  —  2 ax  +  b  y,  sustituyendo  en  la  ecuación,  2x(2ax  +  b)  —3 (ax2  +  bx  +  c)  —  4x2. 
Si  igualamos  los  términos  del  mismo  grado  encontramos  que  esto  se  cumple  con  a  —  4, 
b  —  c  =  0.  Así  pues,  una  solución  particular  de  la  lineal  es  y  =  4a;2.  Por  otra  parte, 
y  tal  como  ya  hemos  calculado  en  (ii),  la  solución  general  de  la  homogénea  asociada 
es  y  =  (7a;3/2.  En  definitiva,  de  nuevo  tenemos  que  la  solución  general  de  la  lineal  es 
y  =  (7  a;3/2  +  4a;2. 

(iv)  Para  resolver  la  ecuación  lineal,  descomponemos  y  =  u(x)v(x),  de  donde  y'  = 
u'v  +  uv' .  Sustituyendo,  2 x(u'v  +  uv')  —  3 uv  —  4a;2,  es  decir,  2 xu'v  +  (2xv'  —  3v)u  =  4a;2. 
Si  igualamos  a  0  el  coeficiente  de  u  queda  2xv  —  3v  —  0,  o  sea  ~  =  integrando, 
log v  —  |loga;,  luego  v  =  a;3/2.  Con  esto,  la  ecuación  queda  2 xu'x3/2  —  4a;2,  es  decir, 
v!  —  2a;-1/2,  cuya  solución  es  u  —  4a;1/2  +  (7.  Sin  más  que  recomponer  y  —  uv  obtenemos 
y  =  (4a;1/2  +  (7)a;3/2,  la  misma  solución  para  la  lineal  que  con  los  otros  tres  métodos. 

APARTADO  5'. 

Ecuación  de  Bernouilli. 

Consideremos 

y’  +  a(x)y  +  b(x)ya  =  0. 

Es  claro  que,  si  a  =  0,  la  E.  D.  anterior  es  lineal  y,  si  a  —  1,  es  de  variables  separadas.  En 
otro  caso,  veamos  cómo  resolverla: 

Efectuemos  el  cambio  de  función  y1-"  =  z,  para  el  cual  (1  —  a)y~ay'  —  z' ,  es  decir, 

—  pry  -  Sustituyendo  en  \y'  +  a(x)y1~°í  +  b{x)  =  0  tenemos  — I -a(x)z  +  b(x)  —  0, 
con  lo  que  hemos  transformado  la  ecuación  de  Bernouilli  en  la  E.  D.  lineal  z'+(l—a)a(x)z  — 
( a  —  l)6(x). 

Puede  seguirse  un  segundo  método  de  resolución  sin  más  que  aplicar  un  procedimiento 
análogo  al  último  de  los  que  hemos  visto  para  ecuaciones  lineales.  Partimos  de  la 
descoposición  y{x)  =  u(x)v(x).  Derivando,  y'  =  u'v  +  uv'  y,  sustituyendo  en  la  E.  D., 
u'v  +  uv'  +  a(x)uv  +  b(x)uava  =  0,  que  escribimos  como  u'v  +  (v'  +  a(x)v)u  +  b(x)uava  —  0. 
Ahora,  igualamos  a  0  el  coeficiente  de  u.  con  lo  cual  tenemos  v'  +  a(x)v  —  0,  que  es  una 
ecuación  en  v(x)  de  variables  separadas;  resolviéndola  determinamos  v(x).  Con  esta  v, 
la  ecuación  de  la  que  partíamos  ha  quedado  u'v  +  b(x)uava  =  0,  que  es  una  ecuación 
en  u(x)  de  variables  separadas.  Resolviéndola  encontramos  u(x),  con  lo  que  ya  tenemos 
completamente  solucionada  la  ecuación  de  Bernouilli. 

Ejemplo  5'.  Resolver 

xy  +  2y  +  x5y3ex  =  0. 

Puesta  en  la  forma  y'  +  ^¡y  +  x4ysex  =  0,  es  claro  que  la  ecuación  es  de  Bernouilli 
con  a  —  3.  Hacemos  el  cambio  de  función  y~ 2  —  z ,  para  el  cual  z'  =  —2y~3y',  es  decir, 
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RECETA  5'.  Ecuación  de  Bernouilli. 

Es  de  la  forma 

y'  +  a(x)y  +  b(x)ya  —  0. 

Si  a  —  0  es  lineal,  y  si  a  —  1,  de  variables  separadas.  En  otro  caso,  se 
hace  el  cambio  de  función  y1-a  =  z,  con  lo  que  la  E.  D.  de  Bernouilli 
se  transforma  en  una  lineal.  Un  segundo  método  de  resolución  es  el 
siguiente:  se  descompone  y(x)  —  u(x)v(x)  y  se  sustituye  en  la  E.  D.,  se 
iguala  a  0  el  coeficiente  de  u  (queda  v'  +  a(x)v  =  0,  que  es  de  variables 
separadas),  lo  que  nos  lleva  a  determinar  v .  apareciendo  ahora  una 
ecuación  en  u(x)  de  variables  separadas. 


^3  =  —  y.  Si  sustituimos  esto  en  ^3  +  f 2/  2  +  x4ex  =  0  queda  z'  —  =  2x4ex,  que  es 

lineal  con  a(x)  —  ~^r  y  b(x)  —  2x4ex. 

Tal  como  sabemos,  la  solución  de  la  ecuación  lineal  es  z  =  exp  (—  f  a(x)  dx)  x 
[/  6(x)exp  (J  a(x)  dx )  dx  +  C]  .  Como  f  a(x )  dx  —  f  ^  dx  —  —4  loga;  =  loga;-4,  se  sigue 
z  —  x4(f  2xxf  dx  +  C)  —  x4(2ex  +  C).  Por  tanto,  la  solución  de  la  E.  D.  de  Bernouilli  es 
y~2  —  x4(2ex  +  C ),  es  decir,  y  —  x~2(2ex  +  C)-1/2 . 

Resolvamos  ahora  la  ecuación  por  el  segundo  procedimiento  explicado.  Tomamos 
y  =  uv,  luego  y'  =  u'v  +  uv'\  sustituyendo,  x(v!v  +  uv')  +  2 uv  +  x5u3v3ex  =  0,  es 
decir,  xu'v  +  (xv'  +  2 v)u  +  x5u3v3ex  =  0.  Elijamos  v  tal  que  xv'  +  2v  —  0,  esto  es, 
1°  cual  se  consigue  con  logn  =:  —2 loga;,  o  sea,  v  =  x~2.  Así,  tenemos  la  nueva 
ecuación  xv!  +  x5u3x~6ex  =  0,  que  puede  escribirse  como  ex  dx.  Su  solución  es 

2^2  —  ex  +  C / 2,  es  decir,  u  =  (2ex  +  C)-1/2.  Con  todo  esto,  la  solución  de  la  ecuación  de 
Bernouilli  de  la  que  partíamos  es  y  —  x~2(2ex  +  C)~1!2.  Lógicamente,  hemos  obtenido  la 
misma  que  mediante  el  primer  método. 


APARTADO  5". 

Ecuación  de  Riccati. 

Supongamos  que  tenemos  la  E.  D. 

y'  +  a(x)y  +  b(x)y2  =  c(x). 

Para  resolverla,  tenemos  que  haber  encontrado  previamente  una  solución  particular  yp(x). 
Si  éste  es  el  caso,  efectuamos  el  cambio  de  función  y  —  yp  +  z,  con  lo  cual  y'  —  y'p  +  z'  y, 
sustituyendo,  yp+z' +a(x)(yp+z)+b(x)(yp+2ypz+z2)  =  c(x).  Como  y'p+a(x)yp+b(x)y2  = 
c(x)  por  ser  yp  solución  particular,  esa  expresión  queda  z'  +  a(x)z  +  b(x)z2  +  2 b(x)zyp  =  0, 


E.  D.  explícitas  de  primer  orden 


23 


es  decir,  z'  +  [a(x)  +  2 b(x)yp(x)\z  +  b(x)z2  —  O,  que  es  una  E.  D.  de  Bernouilli  con  a  —  2. 
(Nótese  ahora  que,  como  el  cambio  z  —  u~x  eu  la  de  Bernouilli  reduciría  ésta  a  una  lineal, 
el  cambio  directo  y  —  yp  +  -  transforma  la  ecuación  de  Riccati  en  una  lineal  de  un  solo 
paso.) 

Existen  diversos  resultados  interesantes  sobre  la  ecuación  de  Riccati.  Por  ejemplo,  si  se 
conocen  dos  soluciones  y\  e  y2  suyas,  el  cambio  de  función  y  =  Vl lleva  a  que  se  pueda 
resolver  la  ecuación  con  una  sola  cuadratura:  v  =  C exp(f  b(x)(y 2(x)  —  yi(x))dx).  Pero 
desafortunadamente,  no  hay  ningún  procedimiento  general  que  permita  obtener  alguna 
solución  particular,  por  lo  que  muchas  veces  las  E.  D.  de  Riccati  no  resultan  integrables 
en  términos  de  funciones  elementales.  A  este  respecto,  resulta  apropiado  citar  la  ecuación 
especial  de  Riccati  y'  +  by2  =  cxm  con  b,  c  G  M  \  {0};  se  sabe  que  puede  resolverse  en 
términos  finitos  si  y  sólo  si  m  —  —  2  o  m  —  para  algún  entero  k. 


RECETA  5".  Ecuación  de  Riccati. 

Es  de  la  forma 

y'  +  a(x)y  +  b(x)y2  —  c(x). 

El  método  requiere  haber  encontrado  previamente  una  solución  parti¬ 
cular  yp(x).  Si  éste  es  el  caso,  haciendo  el  cambio  de  función  y  =  yp  +  z, 
la  E.  D.  de  Riccati  se  reduce  a  una  de  Bernouilli  con  a  =  2. 


Ejemplo  5".  Resolver 


y'  +  y2  —  x2  —  2x. 


Un  primer  vistazo  a  la  ecuación  nos  induce  a  pensar  que,  posiblemente,  tenga  como 
solución  un  polinomio  de  primer  grado  y  —  ax  +  b.  Sustituyendo,  a  +  (ax  +  b)2  —  x2  —  2x, 
lo  que,  efectivamente,  se  satisface  para  a  —  —  1  y  6  =  1.  Por  tanto,  una  solución  particular 
es  yp  =  — x  +  1.  Si  efectuamos  el  cambio  de  variable  y  =  —x  +  1  +  z  (luego  y'  =  —  1  +  z')  y 
sustituimos  en  la  E.  D.,  obtenemos,  tras  simplificar,  z'  +  (2  —  2 x)z  +  z2  =  0,  ecuación  de 
Bernouilli  con  a  =  2.  Haciendo  el  cambio  z  —  u~l  (de  donde  v!  —  Erz'),  al  sustituir  en 
Etz'  +  2~2x  +  1  =  0  aparece  u'  +  (2x  —  2 )u  —  1,  que  es  una  E.  D.  lineal  con  a(x)  —  2x  —  2 
y  b(x)  —  1.  La  solución  de  esta  última  es 


u  =  exp  (  f  (2  —  2x)  dx) 


exp  (/( 2  x  —  2)  dx)  dx  +  C 


=  e 


2x—xz 


^xz  —  2x 


dx  +  C 


Deshaciendo  los  cambios,  u  =  \  ,  luego 


y  +  x  —  1 


—  e 


2x—x2 


^x2  —  2x 


dx  +  C 
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es  la  solución  de  la  E.  D.  de  Riccati.  (De  aquí  puede  despejarse  y  si  se  desea,  obteniéndose 
la  solución  explícitamente,  pero  la  integral  que  aparece  no  puede  resolverse  en  términos  de 
funciones  elementales.) 

APARTADO  6. 

Sustituciones. 

Supongamos  que  tenemos  la  E.  D. 

y'  =  f(x,y) 

que  tiene  un  aspecto  diferente  a  cualquiera  de  las  que  ya  se  lian  estudiado.  En  ocasiones 
puede  ocurrir  que,  sustituyendo  una  parte  de  la  ecuación  por  una  nueva  variable,  la  E.  D. 
aparentemente  difícil  de  la  que  partíamos  se  transforme  en  una  que  puede  resolverse  con 
facilidad.  Esto  es,  en  esencia,  encontrar  un  cambio  de  variables  inteligente.  Aunque  no 
pueden  darse  regias  fijas  sobre  qué  sustituciones  usar,  si  es  que  hay  alguna  sustitución 
posible,  vale  la  pena  intentar  algo  cuando  no  se  nos  ocurre  otro  camino.  Muchas  veces,  una 
sustitución  muy  simple  puede  ser  suficiente. 


RECETA  6.  Sustituciones. 

Cuando  tenemos  una  E.  D. 

y'  =  f(x,y) 

que  no  responde  a  alguno  de  los  tipos  estudiados  hasta  ahora,  a  veces 
una  sustitución  (en  esencia,  un  cambio  de  variable)  más  o  menos 
ingeniosa  transforma  la  ecuación  en  una  reconocible.  Lógicamente,  no 
puede  darse  una  regla  general  pero,  en  todo  caso,  merece  la  pena 
intentar  algo. 


Ejemplo  6.  Resolver 


dy  „  x 3 

x  — - y  —  2 — e 

dx  y 


y/x 


La  ecuación  no  es  separable,  ni  homogénea,  ni  exacta,  ni  admite  factores  integrantes 
fj,(x)  o  n(y),  ni,  finalmente,  es  lineal,  de  Bernouilli  o  de  Riccati.  Sin  embargo,  si  la 
ponemos  en  la  forma  xy  y  —  2-ey^x,  nos  damos  cuenta  rápidamente  que  la  derivada 

x  y 

del  exponente  y/x  es  el  primer  miembro  v  x^v  •  Así,  al  efectuar  el  cambio  de  función 
dado  por  la  sustitución  u  =  -,  la  ecuación  queda  transformada  en  v!  =  2u~1eu,  es 
decir,  ue~u  du  =  2  dx ,  ecuación  en  variables  separadas.  Integrando  por  partes  resulta 
—ue~u  —  e~u  +  C  —  2x,  de  donde  u  +  1  =  (C  —  2x)eu.  Volviendo  a  las  variables  iniciales, 
|  +  1  =  (C  —  2x)ey/x:  y  por  tanto  la  solución  de  la  E.  D.  de  partida  es,  expresada  en  forma 
implícita,  y  +  x  —  x(C  —  2x)ey¡x . 


ECUACIONES  EN  LAS  QUE 

LA  DERIVADA  APARECE  IMPLÍCITAMENTE 


Son  las  que  tienen  la  forma 


F(x,y,y')  =  0. 


APARTADO  7. 

F  algebraica  en  y'  de  grado  n. 

Supongamos  que  tenemos 

(: y')n  +  «i {x:y){y')n~l  H - f  an-1(x,y)y'  +  an{x,y )  =  0. 

Si  pensamos  en  la  expresión  anterior  como  en  un  polinomio  en  y ’  de  grado  n  igualado 
a  cero,  y  logramos  resolver  la  expresión  algebraica,  obtenemos  las  raíces  y'  —  /,(.x,  y), 
i  =  1,  2, . . . ,  n.  Es  decir, 

0 y '  ~  h(x,y))(y'  - h (x,y))---  ( y'  -  fn(x,y ))  =  o. 

Por  lo  tanto,  las  soluciones  de  la  E.  D.  de  partida  serán  las  de  cada  una  de  las  nuevas 
ecuaciones  diferenciales  y'  —  fi(x,y )  =  0,  i  =  1,  2, . . . ,  n,  que  habrá  que  resolver.  De  esta 
forma  obtenemos  n  familias  unipar amétricas  de  soluciones. 


RECETA  7.  F  algebraica  en  y'  de  grado  n. 

Tenemos 

C v'T  +  «i {x,y)(y')n~1  H - f  an-1(x,y)y'  +  an{x,y )  =  0. 

Resolviéndolo  como  un  polinomio  en  y ’  de  grado  n  igualado  a  cero 
obtenemos 

(: y '  -  h(x,y))(,y'  -  h (x,y))---  (y'  -  fn(x,y ))  =  o. 

Por  tanto,  las  soluciones  de  la  E.  D.  de  partida  serán  las  soluciones  de 
cada  una  de  las  ecuaciones  y'  —  fi(x,  y)  =0,  i  =  1,2, ,  n. 
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Ejemplo  7.  Resolver 

y\{y')2  +  1)  =  1. 

Despejando  {y')2  queda  {y')2  =  1~|  ,  cuyas  soluciones  algebraicas  son  y'  =  . 

Vamos  a  resolver  estas  dos  nuevas  ecuaciones  diferenciales  a  la  vez  (ambas  son  de  variables 

separadas).  Si  las  ponemos  como  ±  ^dy  =  dx,  integrando,  =F\/l  —  y2  =  x  +  C.  Si 

v1-y2 

elevamos  al  cuadrado  ambos  términos,  podemos  expresar  conjuntamente  las  dos  familias 
de  soluciones  como  1  — y2  =  (x  +  C)2,  es  decir,  (x  +  C)2  +y2  =  1,  que  son  las  circunferencias 
con  centro  en  el  eje  x  y  radio  1. 

Por  último,  es  evidente  que  y  =  1  e  y  —  —  1  también  son  soluciones  de  la  E.  D., 
aunque  no  se  encuentran  entre  las  que  acabamos  de  hallar.  Si  atendemos  a  la  interpretación 
geométrica  de  las  ecuaciones  diferenciales,  es  lógico  que  estas  dos  rectas  sean  soluciones, 
ya  que  son  las  envolventes  de  la  familia  de  circunferencias  que  satisfacen  la  E.  D. 

Obtención  de  la  envolvente  de  una  familia  de  curvas. 

Recordemos  qué  es  una  envolvente:  Dada  una  familia  uniparamétrica  de  curvas 
Wo]  aeA  en  el  plano,  la  envolvente  de  la  familia  es  una  nueva  curva  ip  tal  que,  en  cada 
punto  de  contacto  de  p  con  alguna  de  las  <pQ,  la  tangente  de  <p  y  de  pa  es  la  misma. 

•  Supongamos  en  primer  lugar  que  tenemos  una  familia  de  curvas  expresadas 
implícitamente  como  F(x,y,C )  =  0  (es  decir,  para  cada  valor  del  parámetro  C  aparece 
una  curva  de  la  familia).  Las  envolventes  se  obtienen  eliminando  C  del  sistema 

í  F(x,y,C)  =  0 
\Fc{x,y,C)=  0  ’ 

donde  con  Fq  hemos  denotado  la  derivada  parcial  de  F  respecto  a  C. 

En  el  caso  de  las  circunferencias  del  ejemplo  anterior  tendríamos 

í  (x  +  C)2  +  y2  =  l 
\  2(x  +  C)  -  0 

Despejando  en  la  segunda  ecuación,  x  =  — C  y,  sustituyendo  en  la  primera,  y2  =  1.  Es 
decir,  las  envolventes  son  las  rectas  y  =  ±1. 

•  En  segundo  lugar,  supongamos  que  tenemos  la  familia  de  curvas  en  paramétricas 

{x  =  x(t,  C ) 

y  =  y(t,C) 

Las  envolventes  se  obtienen  despejando  C  —  C(t )  en  el  jacobiano  igualado  a  cero 

detPc  !,c)=0 

\x,  y¡  J 


E.  D.  en  las  que  la  derivada  aparece  implícitamente 
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y  sustituyendo  en  las  ecuaciones  paramétricas. 

Realmente,  con  estos  procedimientos  no  sólo  aparecen  las  envolventes,  sino  también 
los  lugares  geométricos  de  puntos  singulares ,  ya  sea  de  puntos  de  retroceso,  de  puntos  de 
inflexión  o  de  cúspides.  (En  los  puntos  de  retroceso  no  hay  vector  tangente,  y  surgen  en 
el  proceso  anterior  sea  cual  sea  la  parametrización  de  las  curvas;  en  cambio,  los  puntos 
de  inflexión  y  las  cúspides  pueden  aparecer  o  no  según  sea  la  parametrización  usada.  Por 
ejemplo,  nunca  surgen  al  usar  el  parámetro  arco.  Por  esta  razón,  a  los  primeros  se  les 
llama  puntos  singulares  esenciales  y,  a  los  otros,  puntos  singulares  evitables.)  En  todo 
caso,  después  de  aplicar  estos  métodos,  hay  que  comprobar  siempre  qué  es  lo  que  hemos 
encontrado. 

APARTADO  8. 

Ecuación  de  la  forma  y  =  f(x,y'). 

Como  procedimiento  general  para  intentar  resolver  este  tipo  de  ecuaciones,  tomamos 
y'  —  p  y  derivamos  y  =  f(x,y')  respecto  de  x,  con  lo  cual  tenemos 

P  =  y'  =  fx  +  fy’  ^  =  fx(x,p)  +  fy>{x,p)p'. 

(Quizás  resulte  más  sencillo  interpretar  esto  como  que  derivamos  y  —  f(x. p)  respecto  de 
x,  obteniéndose  directamente  p  —  fx  +  fpp'  ■)  Cuando  /  tiene  la  forma  adecuada,  la  nueva 
ecuación  p  —  fx(x,p)  +  fp(x,p)p'  que  hemos  encontrado  puede  ser  de  alguno  de  los  tipos  ya 
estudiados.  Si  éste  es  el  caso,  la  resolvemos,  obteniendo  su  solución  x  —  (f>(p,  C ).  Entonces, 
la  solución  de  la  E.  D.  de  partida  será 

{x  =  <t>(p,  C) 

\  y  =  f((t>(p,c),p )  ’ 

expresado  como  una  familia  de  curvas  en  paramétricas.  Puede  resultar  extraño  pensar  que 
el  parámetro  p  vale  precisamente  4^;  pero  esto  no  importa  en  absoluto,  sino  que  puede 
considerarse  una  simple  curiosidad. 

Sin  embargo,  no  siempre  ocurre,  para  una  función  /  genérica,  que  el  proceso  anterior 
conduzca  a  una  ecuación  reconocible.  Esto  sí  que  sucede  en  los  tres  casos  que  estudiamos 
a  continuación: 

8.1.  Ecuación 

y  =  f{y')- 

En  este  caso,  como  f(x,y')  —  f(y'),  se  tiene  fx(x,p)  —  0,  luego  con  el  proceso  descrito 
habremos  obtenido  la  ecuación  p  =  fp(p)p'  o,  lo  que  es  lo  mismo,  p  —  f'{p)p' ■  Si  p  ^  0, 

esto  lo  podemos  poner  como  dx  =  ^  dp ,  que  es  de  variables  separadas.  Integrándola, 

x  —  f  ^  dp  —  (f)(p)  +  (7,  y  las  soluciones  de  y  —  f{y')  serán  las  curvas 


x  =  4>(p)  +  C 

y  =  f(p ) 
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Notar,  a  la  vista  de  su  representación  paramétrica,  que  todas  ellas  tienen  la  misma  forma: 
se  diferencian  únicamente  en  un  desplazamiento  horizontal. 

A  partir  de  p  =  0  obtenemos  la  solución  constante  y  —  /(O)  (al  sustituir  en  la  ecuación 
queda  /(O)  =  /(O)).  Gráficamente,  esto  es  una  recta  envolvente  de  las  demás  soluciones. 

8.2.  Ecuación  de  Lagrange: 

V  +  xip(y')  +  =  0. 

Si  ponemos  y'  —  p  y  derivamos  respecto  de  x  queda 

P  +  <p(p)  +  x<p'(p)  ^  +  i/j'ip)  ^  =  0 
dx  dx 

que,  si  lo  preferimos,  podemos  escribir  como  (p  +  <p(p))  dx  +  xip'(p)  dp  +  'i p'{p)  dp  —  0. 

•  Supongamos  en  primer  lugar  que  p  +  <p(p)  ^  0.  Entonces,  podemos  dividir  por 
P  +  (PÍp)i  con  1°  que  nos  aparece  la  ecuación  lineal 

dx  p\p)  x  ip'{p)  =  Q 

dp  p  +  (p(p)  p  +  (p(p) 

en  la  que  x  actúa  como  función  y  p  como  variable.  Si  la  resolvemos,  obtendremos  como 
solución  x  =  (f){p,C).  Con  esto,  habremos  encontrado  para  la  E.  D.  de  Lagrange  las 
soluciones 

(x  =  0(p,  C) 

\y  =  -<¡>{p,c)(pip)  ~^{p) 

en  forma  de  familia  de  curvas  paramétricas. 

•  Supongamos  ahora  que  existe  algún  A  para  el  cual  A  +  ¡p(A)  =  0.  Formalmente, 
tomamos  y1  =  A  y,  sustituyendo  en  la  E.  D.,  y  +  xp( A)  +-¡/>(A)  =  0,  es  decir,  y  =  \x  —  ^(A). 
Esta  recta  es,  efectivamente,  solución  de  la  ecuación  de  Lagrange,  sin  más  que  comprobar 
que  la  satisface:  Aa;  —  ^>(A)  +xtp(\)  +V,(A)  =  0  pues  A  +  <p(A)  =  0.  Estas  rectas  decimos  que 
son  las  soluciones  singulares  de  la  ecuación.  Además,  recíprocamente,  vamos  a  ver  que  si 
y  —  \x  —  es  solución  de  la  E.  D.,  entonces  se  cumple  que  A  +  <p(\)  =  0.  Para  ello,  no 
hay  más  que  utilizar  que  y  =  Xx  —  V’(A)  satisface  la  ecuación,  con  lo  cual  deberá  cumplirse 
0  =  Xx  —  'ip(X)  +  xp(X)  +  ip(X)  —  (A  +  (p(X))x,  luego  A  +  <p(A)  =  0. 

8.3.  Ecuación  de  Clairaut: 


y  -  xy'  +  'ip(y')  =  0. 

Este  es  un  caso  particular  de  ecuación  de  Lagrange  con  <p(y')  —  — y Aquí,  como 
<p( A)  +  A  =  0  para  todo  A  G  R,  no  existen  las  soluciones  que,  en  la  ecuación  de  Lagrange, 
encontrábamos  por  el  método  general;  sólo  aparecen  rectas.  En  realidad,  tenemos  ahora 
como  soluciones  de  la  E.  D.  toda  la  familia  de  rectas  y  —  Xx  —  X ),  A  G  M. 


E.  D.  en  las  que  la  derivada  aparece  implícitamente 
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Vamos  a  ver  que,  además,  existe  una  solución  singular,  la  envolvente  de  este  haz  de 
rectas.  Para  ello,  planteamos  el  sistema 

(  y  =  \X  -  ij)( A) 

1  0  =  x  —  i//(A) 


en  el  que  la  segunda  ecuación  es  la  primera  derivada  respecto  al  parámetro  A.  Para 
encontrar  explícitamente  la  envolvente  habría  que  despejar  A  en  una  de  las  dos  ecuaciones 
y  sustituir  en  la  otra,  lo  cual  no  siempre  es  posible.  Pero  podemos  dejar  su  ecuación 
paramétrica  que,  obviamente,  será 


\x  =  ^'(A) 

\y  =  A^(A)-^(A)  ' 

Aparte  de  la  consideración  geométrica  de  que  las  rectas  no  tienen  puntos  de  retroceso, 
cúspides  o  puntos  de  inflexión,  realmente  podemos  garantizar  que  esta  curva  es  una 
envolvente  y  no  un  lugar  geométrico  de  puntos  singulares.  En  efecto,  en  cada  punto 
(x(X),y(X))  de  la  curva,  esta  interseca  precisamente  a  la  recta  y  —  Xx  —  ip(X)  del  haz; 
y  la  pendiente  de  esa  recta  en  el  punto  de  intersección  (y  en  cualquier  otro,  de  hecho),  que 
es  A,  coincide  con  la  de  la  curva: 

dy  dy/dX  ip' (X)  +  Xip"  (X)  —  ip' (X) 

dx  dx/dX  ^"(X) 

De  paso,  esto  significa  que  las  soluciones  rectas  son  las  tangentes  de  la  solución  singular. 

Realmente,  podemos  hallar  las  soluciones  de  la  ecuación  de  Clairaut  sin  necesidad  de 

saber  que  es  un  caso  particular  de  la  ecuación  de  Lagrange  cuyas  soluciones  son  rectas  y  su 

envolvente.  Para  ello,  tomamos  y'  =  p  e n  la  E.  D.  y  derivamos  la  expresión  y  —  xp+ip(p)  —  0 

respecto  de  x,  obteniendo  p  —  p  —  xp'  +  'i¡r'(p)p'  =  0,  es  decir,  (—x  +  ^'{p))p'  =  0.  De 

,2 

aquí,  p'  —  0  o  ipf(p)  =  x.  En  primer  lugar,  si  suponemos  =  p'  —  0,  tendremos 
^  =  p  —  X  constante;  sustituyendo  en  la  E.  D.  encontramos  la  familia  de  soluciones  rectas 
y  =  Xx  —  ip(X),  AeK.  En  segundo  lugar,  si  suponemos  ^'(p)  —  x,  encontramos  la  solución 
dada  paramétricamente  por 

{x  =  ip\p) 

y  =  p^'íp)  -  VKp) 

(donde,  anecdóticamente,  el  parámetro  p  coincide  con  p  —  ^);  ésta  es  la  solución  singular, 
la  envolvente  del  haz  de  rectas. 

Ejemplo  8.1.  Resolver 

y  =  (y')2  +  2(y')3. 


Tomando  y'  —  p  podemos  poner  y  =  p2  +  2 p3.  Derivando  respecto  de  x  tenemos 
y'  =  p  —  2 pp'  +  6 p2p'  =  (2 p  +  6p2)p' ,  es  decir,  p  —  (2  +  6 p)pp' .  Para  p  ^  0,  simplificamos 
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RECETA  8.  Ecuación  de  la  forma  y  —  f(x,y'). 

En  general,  se  toma  y'  =  p  y  se  deriva  la  ecuación  y  =  f(x,y') 
respecto  de  x.  Si  /  tiene  la  forma  adecuada,  a  la  nueva  E.  D.  se  le 
puede  aplicar  alguno  de  los  métodos  ya  estudiados,  procediéndose  así 
a  su  resolución. 

8.1.  Cuando 

y  =  f(y '), 

la  ecuación  que  se  obtiene  mediante  el  proceso  anterior  es  de  variables 
separadas. 

8.2.  Ecuación  de  Lagrange: 

y  +  x(p(y')  +  ^(y')  =  0. 

Se  reduce  a  una  ecuación  lineal  con  x  como  función  y  p  como  variable. 
Además,  para  los  A  tales  que  A  +  <p(A)  =  0  se  obtienen  como  soluciones 
las  rectas  y  =  Xx  —  V’(A). 

8.3.  Ecuación  de  Clairaut: 

y  -  xy'  +  ijj(y')  =  0. 

Es  un  caso  particular  de  ecuación  de  Lagrange  en  el  que  sólo  aparecen 
rectas  (y  su  envolvente). 


por  p,  con  lo  que  1  =  (2  +  6 p)p' ,  es  decir,  dx  =  (2  +  6 p)  dp.  Integrando,  obtenemos  como 
solución  la  familia  de  curvas  paramétricas 

í  x  —  2 p  +  3 p2  +  C 

\  V  —  P2  +  2 p3 

Además,  a  partir  de  p  —  0  se  tiene  la  solución  singular  y  =  0. 

Para  concluir  el  ejemplo,  vamos  a  ver  qué  obtenemos  si  intentamos  encontrar  la 
envolvente  de  la  familia  de  curvas  paramétricas  anteriores.  Se  tiene 

det(2  +  6p  2 p+V)  =P(2  +  6p)=0  <=>  p  —  0  o  p  — 

Con  p  —  0  obtenemos  la  recta  y  —  0  que,  efectivamente,  es  solución  de  la  E.  D.:  es 
la  envolvente  de  las  otras  soluciones.  Pero  con  p  =  aparece  la  recta  y  —  que, 
claramente,  no  es  solución.  En  realidad,  en  un  punto  (x,y)  de  cualquiera  de  las  curvas,  el 
vector  tangente  es  (x,y)  —  (%,  ^)  =  (2  +  6p,2p  +  6 p2).  Así,  cuando  p  =  0  tenemos  el 
vector  tangente  (i,  y)  =  (2,  0);  pero  cuando  p  —  el  vector  tangente  es  (i,  y)  —  (0,  0),  es 
decir,  no  hay  vector  tangente.  Esto  quiere  decir  que  la  recta  y  —  ^  es  el  lugar  geométrico 
de  los  puntos  de  retroceso  de  la  familia  de  curvas  solución. 


E.  D.  en  las  que  la  derivada  aparece  implícitamente 
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Ejemplo  8.2.  Resolver 

y  =  X  +  y'  -  3 {y')2. 

Nos  encontramos  ante  una  ecuación  de  Lagrange  y+x(p(y')+'i/r(y')  =  0  con  p(y')  —  — 1 
y  %K'yr)  —  ~y¡  +  3{y')2 .  Para  resolverla,  tomamos  p  —  y'  y,  derivando  la  E.  D.  respecto  de 
x ,  tenemos  p  =  1  +p'  —  6 pp',  es  decir,  (p  —  1)  dx  —  (1  —  6 p)  dp.  Para  p  ^  1  podemos  poner 
dx  —  dp  —  (—6—  ^-¡-)  dp,  ecuación  de  variables  separadas.  (Fijarse  que,  al  describir  el 
proceso,  hemos  dicho  que  tenía  que  salir  aquí  una  ecuación  lineal;  por  casualidad,  lo  que  nos 
ha  aparecido  es  una  lineal  homogénea,  que  siempre  es  una  E.  D.  de  variables  separadas.) 
Integrándola,  x  —  —6 p  —  5  log (p  —  1)  +  C.  Sustituyendo  esto  y  y’  —  p  en  la  expresión 
y  =  x  +  y'  —  3(y')2  queda  y  —  —  6p—  5  log(p—  1)  +  C+p  —  3p2  =  —5p  —  3p2  —  5  log(p  —  1)  +  C. 
Así  pues,  las  curvas  paramétricas 

f  x  —  —6 p  —  5  log (p  —  1)  +  C 
|  y  —  —5  p  —  3  p2  —  5  log  (p  —  1)  +  C 

son  soluciones  de  la  ecuación  de  Lagrange  de  la  que  partíamos. 

Además  de  éstas,  a  partir  de  p  —  1  se  obtiene,  sin  más  que  sustituir  literalmente  en 
y  —  x  +  p  —  3 p2,  la  solución  singular  dada  por  la  recta  y  —  x  —  2. 


Ejemplo  8.3.  Resolver 


y  =  y'x-  2 (y')2. 


Estamos  ante  una  ecuación  de  Clairaut  y  —  xy'  +ip(y')  =  0  con  ^{y')  =  2 (y')2.  Hemos 
demostrado  que,  sin  más  que  tomar  y'  =  A,  las  rectas  y  =  \x  —  2A2  son  solución  de  la 
E.  D.  para  cada  A  G  M.  Calculemos  su  envolvente:  en  el  sistema 

í  y  =  Xx  —  2A2 
1  0  =  x  —  4A 


(recordar  que  la  segunda  ecuación  se  obtiene  derivando  la  primera  respecto  al  parámetro  A) 

despejamos  A  =  |  en  la  segunda  ecuación  y,  sustituyendo  en  la  primera,  encontramos  que 

2 

la  envolvente  es  la  parábola  y  = 

Por  otra  parte,  si  hubiéramos  intentado  resolver  la  E.  D.  directamente,  sin  utilizar  lo 
que  hemos  demostrado  previamente,  tendríamos  que  tomar  y '  —  p.  con  lo  cual  y  —  px  —  2 p2 
y,  derivando  respecto  de  x,  p  —  p'x  +  p  —  4 pp' ,  es  decir,  (x  —  4 p)p'  =  0.  De  aquí,  p'  —  0  o 

j2  i 

x  =  4 p.  Si  suponemos  (jc3  =  p'  =  0,  tendremos  ^  —  p  —  A  constante,  luego,  sustituyendo 
en  la  E.  D.,  hubiéramos  encontrado  como  soluciones  las  rectas  y  —  Xx  —  2 A2 ,  A  e  I.  Por 
último,  para  hallar  la  solución  singular  (envolvente  de  las  rectas),  si  suponemos  x  —  4 p, 
tendríamos  la  solución  paramétrica 

{x  —  4  p 

y  —  px  —  2  p2  =  4  p2  —  2  p2  =  2  p2 

lógicamente,  sin  más  que  despejar  p  en  la  primera  expresión  y  sustituir  en  la  segunda,  esta 
curva  es  la  parábola  y  = 


32 


Métodos  clásicos  de  resolución  de  E.  D.  O. 


APARTADO  9. 

Ecuación  de  la  forma  x  =  f(y,y')- 

Estas  ecuaciones  son  similares  a  las  que  aparecen  en  el  apartado  8  pero  con  el  papel 
de  x  e  y  intercambiado.  En  general,  para  intentar  resolverlas,  tomamos  y'  —  p  y  derivamos 
x  =  f(y,y')  (o,  si  se  prefiere  interpretarlo  así,  x  =  f(y,p ))  respecto  de  y  (en  lugar  de 
respecto  a  x,  como  hacíamos  en  el  apartado  8).  Esto  requiere  tener  un  poco  más  de  cuidado, 
además  de  usar  ^  =  ,  ),  —  Así, 

ay  ay/ax  p  7 


i 

p 


dx 

dy 


fy(y,p )  +  fP{y,p ) 


dp 

dy' 


Según  como  sea  /,  esta  nueva  ecuación  ^  =  fy(y,p )  +  fp{y,p )  ^  a  la  que  hemos  reducido 
la  que  teníamos  responde  a  alguno  de  los  tipos  previamente  estudiados,  luego  podríamos 
resolverla.  En  cualquier  caso,  si  logramos  hacerlo,  obtenemos  su  solución  y  =  @(p-  C).  Con 
esto,  la  solución  de  la  ecuación  de  partida  será  la  familia  de  curvas  paramétricas 


x  =  f{<t>(p,C),p ) 
y  =  <t>(p,  c) 


Pero  no  podemos  asegurar  que,  para  una  función  /  cualquiera,  sepamos  resolver  la 
ecuación  intermedia  que  nos  aparece,  con  lo  que  no  habría  forma  de  continuar  el  proceso. 
Se  pueden  estudiar  casos  similares  a  los  de  la  forma  y  —  f(x,y')  en  los  que  se  garantiza 
que  el  método  no  quedará  interrumpido.  Los  tres  casos  que  merece  la  pena  distinguir  son 
los  siguientes: 

9.1.  Ecuación  x  —  f(yr). 

9.2.  Ecuación  x  +  yp{y')  +  il>(yr)  —  0. 

9.3.  Ecuación  x  —  +  ij)(y')  =  0. 

Se  deja  al  lector  que,  como  ejercicio,  describa  el  método  de  resolución  de  estos  tres 
tipos  de  ecuaciones.  Simplemente  hay  que  dedicarse  a  modificar,  con  cuidado,  el  proceso 
utilizado  en  los  correspondientes  tipos  8.1,  8.2  y  8.3,  recordando  que  ahora  hay  que  derivar 
respecto  de  y  en  lugar  de  respecto  de  x. 


RECETA  9.  Ecuación  de  la  forma  x  —  f(y,  y'). 

En  general,  se  toma  y'  —  p  y  se  deriva  la  ecuación  x  —  f(y,y') 
respecto  de  y.  Según  como  sea  /,  la  nueva  E.  D.  que  así  se  obtiene  es 
ya  conocida,  procediéndose  a  su  resolución.  Se  pueden  estudiar  casos 
similares  a  los  de  la  forma  y  =  f(x,y'). 


E.  D.  en  las  que  la  derivada  aparece  implícitamente 
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Ejemplo  9.  Resolver 


x  =  (y'f  +  y'. 


Si  tomamos  =  y1 


p  y  derivamos  x  =  p3  +  p  respecto  de  y  tenemos 


1  _  1 

p  dy/dx 


dx 

dy 


+ 


dp 

dy ’ 


que  podemos  poner  como  dy  =  (3 p3  +p)  dp,  ecuación  en  variables  separadas.  Integrándola, 
y  —  f  PA  +  \p2  +  C.  Por  lo  tanto,  las  soluciones  de  la  ecuación  x  =  (y')3  +y'  son  la  familia 
de  curvas  paramétricas 

{x  —  p3  +  p 

y  =\pA  +  \p2  +  c  ' 


APARTADO  10. 

Ecuación  de  la  forma  F(y,y')  —  0. 

Para  intentar  encontrar  sus  soluciones,  consideremos  la  curva  F(a,/3)  =  0.  El  método 
de  resolución  requiere  que  hayamos  logrado  encontrar  previamente  una  representación 
paramétrica  de  la  curva,  esto  es,  a  =  <p(t)  y  ¡3  —  ip(t)  tal  que  F(<p(t),ip(t))  =  0.  Si  así  lia 
sido,  vamos  ahora  a  efectuar  el  cambio  de  función  y  por  t  mediante  y  —  <p(t),  teniendo  en 
cuenta  que  y'  =  ^(t). 

Entonces,  derivando  y  —  ip(t)  respecto  de  x  tenemos  yr  —  p'{t)  que,  al  ser  y'  =  ip(t), 
puede  escribirse  como 

m  =  ¿(t)  £, 

ecuación  en  variables  separadas. 

•  Si  ip(t)  no  se  anula,  podemos  poner  dx  =  dt  e,  integrando,  x  =  /  dt  +  C. 
Como  consecuencia,  la  familia  de  curvas  paramétricas 

(x  =  Jm-dt  +  c 

1  y  =  <p(t) 

son  soluciones  de  la  E.  D.  de  partida. 

•  Si  existe  algún  í0  para  el  que  —  0,  formalmente  efectuamos  el  siguiente 

razonamiento:  0  =  p'(t)  así  que  p'(t)  =  0  y  consecuentemente  <p(t)  =  cte.  =  p{to), 
lo  que  nos  conduciría  a  la  solución  y  —  ip(to).  Y,  en  efecto,  podemos  comprobar  con 
rigor  que  esta  recta  horizontal  es  solución  de  la  E.  D.  sin  más  que  sustituir  en  ella: 

F{y,  y')  =  F(<p(to),0)  =  F(<p(t0),ip(t0))  =  0. 
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RECETA  10.  Ecuación  de  la  forma  F(y,y')  =  0. 

Consideramos  la  curva  F(a,  ¡ 3 )  =  0.  Si  encontramos  una  represen¬ 
tación  paramétrica  a  —  <p(t),  (5  —  ip{ t ),  F(ip(t),  i/>(t))  =  0,  se  hace  el 
cambio  de  función  y  por  t  mediante  y  —  <p(t),  y'  =  ip(t).  Así,  derivando 
y  —  ip(t )  respecto  de  x.  aparece  una  ecuación  en  variables  separadas. 


Ejemplo  10.  Resolver 

P/s  +  (y')2/3  =  i- 

Si  tomamos  y  —  eos 3t,  y'  =  sen3í,  es  claro  que  (cos3í)2/3  +  (sen3í)2/3  =  1. 
Derivando  y  —  eos3  t  tenemos  y'  =  —  3  eos2  í  sen  í  lo  que,  usando  y'  —  sen3 1.  resulta  ser 
sen3 1  —  —  3  eos2  t  sen  t  ^  ■ 

ax 

Si  suponemos,  en  primer  lugar,  senf  ^  0,  tenemos  dx  —  ~3 1  dt  —  — 3cot 2tdt. 
Integrando  x  =  —3  J  cot2  tdt  —  —3  f  (1  +  cot2  t)  dt  +  3  J  dt  —  3  cot  t  +  3t  +  C.  Por  tanto, 
las  curvas  paramétricas 

(  x  —  3  cot  t  +  3t  +  C 
1  y  —  eos3 1 

son  soluciones  de  la  E.  D.  de  partida. 

Por  último,  para  los  t  tales  que  sen3  t  —  0.  lo  cual  ocurre  cuando  t  =  kx.  k  G  Z, 
aparecen  como  soluciones  singulares,  tal  como  justificamos  al  desarrollar  el  método,  las 
rectas  horizontales  y  =  cos3(&7r),  es  decir,  y  =  1  e  y  =  —1. 


ECUACIONES  DIFERENCIALES  EN  LAS  QUE 
SE  PUEDE  REDUCIR  EL  ORDEN 


Supongamos  que  tenemos  la  E.  D. 

F(x,  y,  y', . . . ,  y^)  =  0  con  n>  1. 

Al  contrario  de  lo  que  ocurría  con  ecuaciones  de  primer  orden,  no  hay  muchos  métodos 
para  encontrar  soluciones  de  ecuaciones  diferenciales  de  orden  mayor  que  uno.  Por  ejemplo, 
está  demostrado  que  la  E.  D.  lineal  general  de  orden  n  >  1  no  es  soluble  por  cuadraturas. 
Sí  que  pueden  resolverse  las  ecuaciones  lineales  de  orden  n  con  coeficientes  constantes; 
aunque  éste  es  un  tema  central  en  el  estudio  de  la  teoría  de  E.  D.  O.,  no  nos  ocuparemos 
de  ello  en  estas  notas. 

Realmente,  lo  único  que  vamos  a  ver  aquí  es  una  serie  de  métodos  que  permiten  reducir 
el  orden  de  una  E.  D.  Aplicándolos  (sucesivamente  si  es  necesario),  podremos  llegar  a  una 
ecuación  de  primer  orden  que,  con  un  poco  de  suerte,  se  encontrará  entre  las  que  ya  sabemos 
resolver. 


APARTADO  11. 


Ecuación  de  la  forma  F(x:y^k\  . . .  ,  y(n))  =  0  con  1  <  k  <  n. 

Tenemos  una  ecuación  en  la  que  no  aparece  la  variable  dependiente  y  (y  además,  si 
k  >  1,  tampoco  sus  derivadas  hasta  el  orden  k  —  1).  Es  evidente  que  el  cambio  y(k^  —  z  la 
convierte  en 

F(x,z,...,zín-k>)  =  0, 

que  es  una  E.  D.  de  orden  n  —  k.  Si  logramos  resolverla,  obtenemos  que  su  solución  será 
z  —  4>(x,  Ci, . . . ,  Cn-k),  una  familia  dependiente  de  n  —  k  constantes.  Entonces,  como 
y(k)  _  para  encontrar  las  soluciones  de  la  ecuación  original  bastará  con  integrar  veces 
( t>(x ,  C i,  •  •  • ,  Cn-k )  respecto  de  x,  es  decir, 


y  = 


(f>(x,  C\, . . . ,  Cn-k )  dx  dx  . . .  dx  . 


k  veces 


k  veces 


lo  cual  introducirá  k  nuevas  constantes  en  la  solución  general.  De  este  modo,  dicha  solución 
general  tendrá  la  forma 


y=  //•/  Cu  •  •  •  5  Cn—k)  dx . . .  dx^  +  Cn-k- \~\X  +  •  •  •  +  Cn—\X  +  Cn. 

k  veces 


k  veces 
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RECETA  11.  Ecuación  de  la  forma  F(x,y(k\  . . .  ,  =  0. 

Mediante  el  cambio  y (fc)  =  z  se  convierte  en  una  ecuación  de 
orden  n  —  k. 


Ejemplo  11.  Resolver 


//  ///  ,  /  ni\ 3  n 

V  -xy  +{y  )  =  0. 


Tomemos  y"  —  z ,  con  lo  que  la  ecuación  se  transforma  en  z  —  xz'  +  (V)3  =  0,  que  es 
de  Clairaut.  Sin  más  que  sustituir  z'  —  A,  las  soluciones  de  ésta  son  las  rectas  z  —  Xx  —  A3, 
AeK.  Además,  también  son  solución  las  envolventes.  Para  calcularlas,  en  el  sistema 

z  —  Xx  —  X3 
0  =  x  -  3X2 


despejamos  A  =  iy'l  en  la  segunda  ecuación,  con  lo  que,  sustituyendo  en  la  primera, 
encontramos  las  dos  envolventes  z  —  x/f'-  Una  vez  resuelta  completamente  la  ecuación 

de  Clairaut,  hallemos  las  soluciones  de  la  original: 

Como  y"  =  z,  a  partir  de  z  —  Xx  —  A3,  integrando  dos  veces, 


y  = 


( Xx  —  A3)  dx  )  dx  —  — - — — f  C\x  +  C2. 

6  2 


Cambiando  la  notación  de  las  constantes,  y  —  K i  x3  —  108 K3x2  +  K¿x  +  K-¿. 
Por  último,  de  z  =  ±%-\/¥  obtenemos 


3  V  3 

±2  3/2  ,  \  ,  f  (  ±2  X 5/2  A  ±4  X7/2 

y  =  I  |  /  zt~híx  dx  )  dx  =  /  (  +  Ci  )  dx  =  +  ^2- 


3^/3" 


V3V/3  5/2 


15V3  7/2 


Es  decir,  las  dos  familias  de  curvas  y  =  +  Cix  +  C2  y  y  —  —  fíyr\/¥  +  Cix  +  C2. 


APARTADO  11'. 

Ecuaciones  lineales  de  orden  superior. 

Supongamos  que  tenemos  una  ecuación 


an(x)y (n)  +  an-i(x)y('n  1}  H - h  ai(x)y/  +  a0(x)y  =  fif(x) 


con  n  >  1.  En  primer  lugar,  conviene  dejar  constancia  de  que  existe  una  teoría  general 
bien  establecida  sobre  ecuaciones  diferenciales  lineales  de  cualquier  orden.  Pero  éste  es 
un  estudio  esencialmente  teórico:  110  es  posible  describir  la  solución  general  de  una  E.  D. 


E.  D.  en  las  que  se  puede  reducir  el  orden 
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lineal  de  orden  2  o  superior  por  medio  de  cuadraturas.  Sin  embargo,  sí  que  se  calculan  las 
soluciones  de  forma  muy  satisfactoria  cuando  nos  encontramos  ante  ecuaciones  lineales  con 
coeficientes  constantes  (es  decir,  cuando  a¿(x)  =  a¿  G  R,  0  <  i  <  n).  Pero  todo  esto  queda 
fuera  de  los  objetivos  que  estamos  persiguiendo  aquí.  Fundamentalmente,  en  este  apartado 
nos  preocuparemos  sólo  de  ver  cómo  una  ecuación  lineal  de  orden  n  >  1  puede  reducirse  de 
orden,  sucesivamente  si  es  necesario,  hasta  llegar  a  una  ecuación  lineal  de  primer  orden,  de 
las  que  liemos  estudiado  en  el  apartado  5.  En  realidad,  este  método  de  reducción  de  orden 
radica  en  haber  encontrado  previamente  una  solución  particular  de  la  lineal  homogénea 
asociada 

an(x)y(n)  +  an-1(x)y(n-1)  H - b  a1(x)y'  +  a0(x)y  =  0; 

si  no  hemos  sabido  hallar  tal  solución  particular,  el  proceso  aquí  descrito  no  tiene  para 
nosotros  ninguna  utilidad  práctica. 

Supongamos  pues  que  hemos  logrado  encontrar  una  solución  particular  yn(x)  de  la 
homogénea  asociada  (más  adelante  se  pondrá  de  manifiesto  por  qué  resulta  adecuado 
introducir  una  n  en  la  notación  de  la  solución  particular).  En  este  caso,  efectuamos  el 
cambio  de  función  y  por  z  dado  por  y{x)  =  yn(x)z(x).  Sin  más  que  derivar  sucesivamente, 
obtenemos 


y'{x)  =  y'n(x)z(x)  +  yn(x)z'(x) 

y"(x)  =  y'ri  (x)z(x)  +  2  y'n(x)z'(x)  +  yn(x)z"(x) 


y{n\x)  =  vt\x)z(x)  +  ■  ■  ■  +  í"j9i‘)(x)z'"-‘>(1)  +  ■  ■  ■  +  yn(x)zin\x). 

Sustituyendo  en  la  lineal,  y  sacando  factor  común  los  z^k\  queda 

anynZ{n)  +  (■ nany'n  +  an_iyn)^(n_1)  H - b  {nany^~l)  -\ - b  2 ci2y'n  +  atyn)^ 

+  (an2/in)  H - b  ai y'n  +  a0yn)z  =  g{x). 

El  coeficiente  de  z  en  la  expresión  anterior  es  0  por  ser  precisamente  yn  solución  de  la 
ecuación  lineal  homogénea.  Así,  podemos  poner 

bn(x)z^n?>  +  6n_i(x)z(n_1)  H - b  bi(x)z '  =  g(x). 

Esta  es  una  E.  D.  de  las  estudiadas  en  el  apartado  anterior.  Si  en  ella  hacemos  z'  —  u 
aparece 

6n(x)'U(n_1)  +  bn-i(x)u('n~ 2-)  H - b  bi(x)u  =  g(x), 

que  es  una  ecuación  lineal  de  orden  n  —  1. 

Si  sabemos  resolver  esta  nueva  ecuación  (por  ejemplo,  si  n  —  1  =  1  y  aplicamos 
cualquiera  de  los  métodos  de  resolución  de  E.  D.  lineales  de  primer  orden  ya  estudiados), 
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tendremos  que  su  solución  general  será  una  función  u(x)  —  (f)(x,  . . . ,  Cn- 1)  dependiente 

de  n  —  1  constantes.  Entonces,  la  solución  de  la  de  partida  será 


y(x)  =  yn(x)z(x)  =  yn(x)  /  u(x)dx  =  yn(x )  /  <j>(x,  Cu  . . . ,  Cn-i)  dx 


=  yn(x )  /  (p(x,Ci, . . .  ,Cn-i)  dx  +  Cnyn(x), 


donde  Cn  surge  como  constante  de  integración  al  calcular  una  primitiva  de  ó  (que  aparece 
representada  con  la  misma  notación  de  integral,  esperando  que  este  pequeño  y  habitual 
abuso  de  notación  no  introduzca  confusión  en  el  lector).  Si  no  conseguimos  resolverla  en 
general,  se  puede  intentar  reducir  el  orden  de  nuevo  encontrando  ahora  un-i(x)  solución 
particular  de  la  homogénea  asociada.  Si  este  proceso  lo  logramos  hacer  el  suficiente  número 
de  veces,  llegaremos  siempre  a  una  ecuación  lineal  de  orden  1,  que  sí  que  sabemos  resolver. 

Finalizaremos  comentando  cómo  varios  de  los  resultados  estudiados  para  E.  D.  lineales 
de  primer  orden  también  pueden  aplicarse  a  lineales  de  orden  superior.  En  primer  lugar,  es 
inmediato  comprobar  (tal  como  hacíamos  en  (m)  del  apartado  5)  que  la  solución  general 
de  la  lineal  puede  expresarse  como  una  solución  particular  de  la  lineal  más  la  solución 
general  de  la  lineal  homogénea  asociada.  (Lógicamente,  si  la  lineal  ya  es  homogénea,  puede 
tomarse  como  solución  particular  la  función  nula.) 

Por  otra  parte,  el  proceso  de  reducción  de  orden  que  hemos  utilizado  nos  muestra,  por 
inducción,  que  la  solución  general  de  la  lineal  de  orden  n  puede  expresarse  como 

y(x )  =  yp(x)  +  Ciyi(x)  H - b  Cn-iyn-í(x)  +  Cnyn(x), 

donde  yp(x)  es  una  solución  particular  de  la  lineal.  En  efecto,  para  orden  n  =  1  lo  hemos 
probado  al  tratar  las  lineales  de  primer  orden.  Y,  para  efectuar  el  paso  de  inducción  de 
orden  n  —  1  a  n,  basta  fijarse  en  el  desarrollo  del  proceso  de  reducción  seguido.  En  efecto, 
por  hipótesis  de  inducción,  la  solución  de  la  E.  D.  en  u(x)  de  orden  n—1  que  nos  aparecía 
tendrá  la  forma 

u(x)  =  <f>(x,  Ci,  .  .  .  jCn-i)  =  Up(x)  +C1U1{x)  H - f  Cn-xUn-xix). 

Por  tanto,  la  solución  de  la  ecuación  de  partida  de  orden  n  será 


y(x)  =  yn(x)z(x)  =  yn(x)  /  u(x)  dx  =  yn{x)  /  C1, . . . ,  Cn_i)  dx 


=  Vn(x)  /  {up{x)  +  CxUlix)  d - f  Cn-iUn-xix))  dx 


=  yn{x)  (  /  up(x)  dx  +  C\  /  ui(x)  dx  +  ■  ■  ■  +  Cn-i  /  un-i(x)  dx  +  CT 


yn{x){Up{x)  CiUi(x')  “b  ‘  ‘  “b  Cn—iUn—\{x)  -\-  Cn) 
=  yP{x)  +  Cxy^x)  H - b  Cn-iyn-x^x)  +  Cnyn(x), 


E.  D.  en  las  que  se  puede  reducir  el  orden 
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tal  como  queríamos  comprobar. 

También  existe  el  método  de  variación  de  las  constantes  para  E.  D.  lineales  de  orden 
n,  y  está  destinado  a  resolver  la  ecuación  lineal  cuando  se  conoce  la  solución  general  de  la 
homogénea  asociada,  y(x )  =  C\yi(x )  +  C2y2(x)  +  •  •  •  +  Cnyn(x).  En  estas  condiciones,  se 
busca  la  solución  de  la  lineal  de  la  forma 

y(x)  =  C1{x)y1{x)  +  C2(x)y2(x)  H - b  Cn(x)yn(x) 

para  ciertas  funciones  Ck{x)  a  determinar.  Para  encontrar  estas  C4(.x),  se  deriva 
sucesivamente  la  expresión  anterior  y  se  sustituye  en  la  ecuación  lineal,  igualando  a  cero 
los  coeficientes  de  los  yk(x).  Aunque  no  nos  preocuparemos  de  ello,  puede  demostrarse 
rigurosamente  que  esto  siempre  conduce  a  la  solución  buscada. 


RECETA  11'.  Ecuaciones  lineales  de  orden  superior. 

Son  de  la  forma 

an{x)y^n)  +  an_i(»y(ri-1)  H - b  ai(x)y'  +  a0(x)y  =  g(x). 

Si  logramos  encontrar  alguna  solución  yn(x)  de  la  lineal  homogénea 
asociada,  el  cambio  de  función  y  —  ynz  hace  que  la  lineal  se  transforme 
en  una  del  tipo  anterior,  cuyo  orden  se  puede  reducir.  Así,  aparece  una 
nueva  ecuación  lineal,  esta  vez  de  orden  n  —  1. 


Ejemplo  11'.  Resolver 

xy"  +  ( 7x  —  1  )y'  —  7  y  —  x2e~7x 

sabiendo  que  la  lineal  homogénea  asociada  tiene  una  solución  de  la  forma  eax . 

En  primer  lugar,  determinemos  a  para  que  y2  ~  eax  sea  solución  de  xy"  +  ( 7x  —  1  )y'  — 
7 y  =  0.  Derivando,  y2  —  aeax ,  y2  —  a2eax;  sustituyendo,  ci2xeax  +  (7x  —  1  )aeax  —  7eax  =  0, 
es  decir,  [x(a2  +  7a)  —  a  —  7}eax  =  0,  lo  cual  se  consigue  con  a  =  —7. 

Apliquemos  ahora  el  método  de  reducción  de  orden.  Tomamos  y  —  y2z,  con  lo  cual, 
derivando,  y'  =  y2z  +  y2z' ,  y"  =  y2z  +  2y2z'  +  y2z" .  Si  sustituimos  en  la  lineal,  queda 
(y"z  +  2y'2z'  +  y2z")x  +  (7x  —  1  )(y2z  +  y2z')  —  7y2z  —  x2e~7x .  Sacando  factor  común  z",  z ' 
y  z,  aparece  xy2z"  +  [2xy2  +  (7x  —  l)y2\z'  +  [xy"  +  (7x  —  l)y2  —  7y2]z  —  x2e~7x.  Como  el 
coeficiente  de  z  coincide  con  la  lineal  homogénea,  de  la  que  y2  es  solución,  esta  ecuación, 
tras  sustituir  y2  e  y2,  se  transforma  en  xe~7xz"  +  \—i4xe~7x  +  (7x  —  1  )e~7x]z'  —  x2e~7x. 
Simplificando  e~7x ,  podemos  poner  xz"  +  (— 7x  —  1  )z'  =  x2. 

(Nótese  que  suele  ser  más  cómodo  utilizar  el  símbolo  y2  en  lugar  de  su  expresión 
concreta  ya  que,  de  este  modo,  no  hace  falta  arrastrar  tampoco  las  expresiones  de  y2  ni 
y 2 ,  sino  que  basta  con  dar  su  valor  al  final  de  las  simplificaciones.  En  particular,  no  se 
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necesita  calcular  yí¿.  Por  el  contrario,  si  empleamos  directamente  las  expresiones  de  y2,  y'2 
e  y'2  desde  el  principio,  esto  puede  ayudarnos  a  detectar  algún  error  que  hayamos  podido 
cometer:  es  fácil  que  el  error  nos  llevara  a  un  punto  sin  salida  en  el  que  110  se  podría 
continuar  el  proceso  al  no  anularse  el  coeficiente  de  z.) 

Tenemos  ahora  una  E.  D.  de  segundo  orden  en  z(x)  en  la  que  no  aparece  explícitamente 
z.  Para  reducirla  de  orden,  basta  tomar  z'  —  u,  con  lo  cual  queda  xu'  +  (— 7x  —  l)it  =  x 2, 
esto  es,  u'  +  (— 7—  -)u  =  x,  lineal  de  primer  orden.  Resolvámosla  aplicando  directamente  la 
fórmula  general  que  hemos  deducido  en  el  apartado  5.  Así,  sin  más  que  tomar  a(x)  —  —7—  - 
y  b(x)  —  x,  calculando  previamente 


exp 


a(x)  dx  )  =  exp 


exp(— 7x  —  loga;)  =  x  le  7x 


tendremos  que  su  solución  es 


—  1—7* 


dx  +  C\  =  xe7x  (-f^e  7x  +  C\)  —  -j^x  +  C\xe‘ 


Deshaciendo  los  cambios, 

y(x)  =  y2(x)z(x)  =  y2(x)  j  u(x)  dx  -  e”7*  J  (=±x  +  Cxxe7x)  dx 

—  e~7x  +  C\  J  xe7x  dx^j  —  e~7x  +  f C\xe7x  —  C±  J  \ e7x  dx^j 

=  e~7x  {^x2  +  ±ClXe7x  -  ±Cie7x  +  C2)  -  ^x2e~7x  +  C,  (±x  -  ¿)  +  C'se-7*. 

Cambiando  la  notación  de  C\  por  49(7i,  la  solución  general  de  la  lineal  de  segundo  orden 
queda 

y{x)  =  =jx2e~7x  +  Ci(7x  -  1)  +  C2e~7x . 


APARTADO  12. 

Ecuación  de  la  forma  F(y,  y' , . . . ,  y =  0. 

En  esta  ecuación  no  aparece  explícitamente  la  variable  independiente  x.  Vamos  a  ver 
cómo  se  reduce  de  orden  si  efectuamos  el  cambio  y’  —  p  y  la  transformamos  en  una  nueva 
en  la  que  aparecerán  y,  p  y  las  derivadas  de  p  respecto  de  y.  (¡Atención!;  respecto  de  y , 
no  respecto  de  x.)  Para  ello,  vamos  a  ver  el  proceso  a  seguir  para  escribir  las  derivadas 
sucesivas  de  y  respecto  de  x  en  función  de  y,  p,  y  las  derivadas  de  p  respecto  de  y.  Ya 
tenemos  —  y'  =  P-  Para  las  siguientes,  si  denotamos  p'  —  p"  =  . . . ,  y  hacemos 

uso  de  la  regla  de  la  cadena,  entonces 

n  dy’  dp  dp  dy  , 

V  =  ~ch¡  =  dx  =  ~dyáx=PP: 


E.  D.  en  las  que  se  puede  reducir  el  orden 
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y 


/// 


dx 


dy_ 

dx 


[p"p  +  {p')2]  p, 


y  así  sucesivamente.  Nótese  que  en  la  expresión  de  cada  y ^  sólo  aparecen  derivadas  de  p 
hasta  el  orden  k  —  1. 

Con  esto,  sustituyendo  en  F  los  valores  que  hemos  encontrado  para  y' ,  y" ,  . . . ,  y (n\ 
la  E.  D.  original  se  transforma  en  una  de  la  forma 


G 


í  dp  dn  1p  \ 
V  ’  dy1 "  '  ’  dyn~l) 


0, 


cuyo  orden  es  n  —  1. 

Si  logramos  resolver  esta  nueva  ecuación  encontraremos  que,  en  general,  su  solución 
será  de  la  forma  p  —  rf)(y,  Ci, . . . ,  (7n_i),  dependiente  de  n  —  1  constantes.  Ahora, 
devolviendo  a  p  su  valor  original,  es  decir,  p  =  obtenemos  jé  =  (p(y,  C i, ,  Cn- 1),  que 
es  otra  ecuación,  esta  vez  de  primer  orden,  que  también  hay  que  resolver  (lo  cual,  además, 
añade  una  nueva  constante).  Las  soluciones  de  esta  última  serán,  obviamente,  las  mismas 
que  las  de  la  E.  D.  de  partida. 


RECETA  12.  Ecuación  de  la  forma  F(y,y,1 . . .  =  0. 

Hacemos  el  cambio  y'  —  p  y  transformamos  la  E.  D.  en  una  nueva 
dependiendo  de  y,  p  y  las  derivadas  de  p  respecto  de  y.  Esta  es  de 
orden  n  —  1. 


Ejemplo  12.  Resolver 


2  yy"  =  0 v T  + 1- 


Tomando  y'  =  p,  y  tal  como  liemos  comprobado  al  explicar  el  método,  se  tiene  y"  =  p'p 
donde  p'  =  Entonces,  al  sustituir  en  la  E.  D.,  ésta  se  nos  transforma  en  2 yp'p  =  p2  +  1, 

es  decir,  que  es  de  variables  separadas.  Si  integramos,  log(p2  +  1)  =  log(C'iy), 

de  donde  p2  +  1  =  C\  y.  y  por  tanto  p  —  ±y/C\y  —  1.  Volvemos  ahora  a  restaurar  el  valor 
p  =  con  lo  cual  nos  aparece  -A  =  ±\/Ciy  —  1.  Esta  es  de  nuevo  una  ecuación  en 
variables  separadas,  que  podemos  escribir  como  (C\y  —  1 ) ~ 1  2  dy  —  ±dx.  Integrándola, 
2^ClV  -  1  =  ±C tx  +  C2,  con  lo  que  ya  hemos  encontrado  las  soluciones  de  la  E.  D.  de 
partida. 
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APARTADO  12'. 

Si  la  ecuación  F(x,  y,  y' , . . . ,  y —  0 
es  tal  que,  para  a  y  m  fijos,  F  cumple 

F(  \x,  Xmu0,  Xm~1ui, ...,  X  m~nun)  =  XaF(x,  u0,u  i, . . . ,  un), 


vamos  a  ver  que,  efectuando  un  cambio  tanto  de  variable  independiente  como  de 
dependiente,  esta  E.  D.  se  podrá  transformar  en  una  del  tipo  anterior  (es  decir,  en  la 
que  no  aparecerá  la  variable  independiente).  Una  ecuación  de  estas  características  suele 
decirse  que  es  homogénea  generalizada  de  grado  a,  en  la  que  cada  factor  x  contribuye  con 
grado  1,  cada  y  con  grado  m,  y '  con  grado  m  —  1,  y"  con  m  —  2,  etcétera. 

Tomamos  dos  nuevas  variables  t  y  z  (t  la  independiente  y  2  la  dependiente)  que 
relacionamos  con  las  originales  x  e  y  mediante  x  —  et ,  y  —  emíz.  Con  esto,  vamos  a  ver 
cómo  escribir  y1  —  ,  y ^  en  función  de  t,  z  y  las  derivadas  de  z  respecto  de  t. 

Si  denotamos  z'  —  4f  (y,  por  supuesto,  z ^  =  4rf ),  tenemos 


áy_  _  dy/dt  _  rnemtz  +  emtz'  _  , 

dx  dx/dt  e1- 


A  partir  de  aquí, 
d2y 


dx2 


d  f  dy 
dx  V  dx 


d  ( dy\  dt  d 


7  \  7  /  7  \e("m  1^í( z'  +  mz )) 

dt  \dx  J  dx  dt  \  / 

(■ m  —  1  )e('m“1-)í(^/  +  mz)  +  e<<rn~l^t{z"  +  mz') 

[z"  +  (2 m  —  1  )z'  +  m(m  —  1  )z) ; 


dx  /  dt 


—  e(m-2  )t(vn 


análogamente 


d3y  d  í  d2y\  d  í  d2y\  dt 
dx 3  dx  \d2x J  dt  \d2x J  dx 

=  (■ m  —  2  )e('rn~2^t(z''  +  (2  m  —  1  )z'  +  m(m  —  1  )z) 

+  e("*-2)*(z"'  +  (2 m  -  1  )z"  +  m(m  -  l)z')]  ^ 

J  e* 

=  eí"1-3)*  (z"'  +  (3  m  —  3  )z"  +  (3m2  —  6  m  +  2  )z'  +  m(m  —  1  ){m  —  2  )z), 


y  así  sucesivamente.  En  general,  por  inducción,  es  claro  que 


d^y 

dx(kl 


z,  z 


Ák) ) _ 


E.  D.  en  las  que  se  puede  reducir  el  orden 
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De  este  modo,  sustituyendo  en  la  E.  D.  que  estamos  intentando  resolver,  obtenemos 
F(é,  emtz,  e^-1)  V  +  mz), ...,  ¿m~n)tgn(z,  z',...,  z^))  =  0; 
extrayendo  X  —  e1,  esto  resulta  ser 

eQÍF(l,  2,  (z '  +  mz), . . . ,  gn(z,  z! , . . . ,  z(n)))  =  0. 


Como  eat  no  puede  anularse,  el  otro  factor  tiene  que  ser  cero,  luego  liemos  transformado 
la  ecuación  de  partida  en  una  de  la  forma 

G(z,z',...,z^)  =  0, 


en  la  que  no  aparece  explícitamente  la  variable  independiente  t.  Esta  ecuación  puede 
reducirse  de  orden  aplicando  el  proceso  descrito  en  el  apartado  anterior.  En  cualquier  caso, 
si  logramos  resolverla,  su  solución  será  z  =  (f>(t,  C\, . . . ,  Cn);  deshaciendo  los  cambios  x  —  e1 
e  y  —  emtz ,  tendremos  que  la  solución  de  la  E.  D.  original  será  y  =  xm(j)( loga;,  C\, . . . ,  Cn). 


RECETA  12'.  Si  la  ecuación  F(x,  y,  y', . . . ,  y^)  —  0 

es  tal  que,  para  a  y  m  fijos,  F  cumple 

F(Xx,  Xmu0,  Xm~1ui, Xm~nun)  =  XaF(x,  u0,  uu  . . . ,  un), 

haciendo  el  cambio  x  =  eL,  y  =  emtz  la  E.  D.  se  transforma  en  una 
de  la  forma  G(z,  z' , . . . ,  z^)  —  0,  a  la  que  se  puede  aplicar  el  método 
anterior. 


Ejemplo  12'.  Resolver 


A  2  3  // 

Ax  y  y 


Analicemos  en  primer  lugar  cómo  lograr  que  esta  ecuación  sea  homogénea  generalizada 
en  la  que  cada  factor  x  contribuya  con  grado  1,  y  los  factores  y,  y',  y"  con  grados 
m,  m  —  1  y  m  —  2  respectivamente.  Para  ello,  cada  monomio  tiene  que  ser  del  mismo 
grado.  En  el  miembro  de  la  derecha,  esto  se  consigue  con  2  =  4m,  es  decir,  m  = 
con  lo  cual  x2  —  y4  es  de  grado  2.  Entonces,  el  miembro  de  la  izquierda  es  de  grado 
2  +  3 m  +  (m  —  2)  =  2  +  §  +  ( |  —  2)  =  2,  coincidente  con  el  del  otro  miembro.  Así  pues,  con 
m  =  \  la  E.  D.  es  homogénea  generalizada  (de  grado  2,  aunque  esto  no  tiene  importancia). 


Hagamos  ahora  los  cambios  x  =  e¿,  y  —  él2z,  y  denotemos  z'  —  4f ,  z"  —  4rf .  De  este 
modo, 

,  dy  dy/dt  \et!2z  +  é^z'  ,/2  , 

=  i  =  £jit  = 1 — ? - = e  ’F  +  b). 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


Sustituyendo  en  la  E.  D.,  Ae2te3tt2z3e~3tl2(z"  —  \z)  =  e2t  —  e2tz4,  de  donde,  simplificando, 
4 z3z"  —  1,  ecuación  en  la  que  uo  aparece  la  variable  independiente  t.  Para  resolverla, 
tomamos  z'  =  p,  con  lo  que 


n  dz  dp  dp  dz  , 

z  —  — —  =  —  =  —  —  —  p  p 

dt  dt  dz  dt 

con  la  notación  p'  —  j^.  Así,  tenemos  la  nueva  ecuación  4 z3p'p  =  1,  con  variable 
dependiente  p  e  independiente  z,  que  ya  es  de  primer  orden.  Esta  ecuación,  que  podemos 
poner  como  4 pdp  =  z~3  dz,  es  de  variables  separadas.  Integrando,  2 p2  =  —  t,z~2  -hC¡ ,  o  sea, 
4p2  =  2Ci  —  z~2 .  Como  p  —  z'  —  %,  hemos  llegado  ahora  a  plantear  4(V)2  =  2C\  —  z~2, 
es  decir,  2 z'  =  ±\J2C\  —  z~2 ,  que  podemos  escribir  como  ±2(2Ci  —  z-2)-1/2  dz  =  dt.  De 
nuevo  integrando, 


t  +  C2  =  ±  I  2(2Ci  -  z~2)~1/2  dz  —  ±  1 2z(2C\Zi  -  1)' 


Wdz 


1  {2C\z2  -  l)1/2  ,  1 


2C*i 


1/2 


=  ±v=r(2C'1^2  —  l)1/2. 


Volviendo  a  las  variables  iniciales  x  —  el  e  y  =  e¿/2^,  aparece  ahora  loga;  +  (72  = 
±^r-(2(7ip2a;-1  —  l)1^2,  con  lo  cual,  sin  más  que  cambiar  la  notación  de  las  constantes, 

K i  loga;  +  K‘2  —  í\/2K ¡ y2x~ 1  —  1 .  Así  pues,  elevando  al  cuadrado  y  despejando  y2,  las 
soluciones  de  la  E.  D.  de  partida  vienen  dadas  por 

2  (Ki  loga;  +  K2)2  +  1 
V  = - 2 K¡x - ’ 

donde  K\  y  K2  son  constantes  arbitrarias. 


APARTADO  13. 

Si  la  ecuación  F(x,  y,y' , ,  y^)  =  0 
es  tal  que,  para  a  fijo,  F  cumple 

F(x,  Xu0,  Xui, . . . ,  A  un)  =  XaF(x,  u0,ui,...,  un), 

vamos  a  comprobar  cómo,  a  través  de  un  cambio  de  variable  dependiente,  el  orden  se 
puede  reducir  en  uno.  La  propiedad  que  caracteriza  a  la  función  F  suele  expresarse 
diciendo  que  F  es  homogénea  de  grado  a  respecto  a  los  argumentos  uq,  ui,  . . . ,  un.  Con  este 


E.  D.  en  las  que  se  puede  reducir  el  orden 
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lenguaje,  la  ecuación  que  estamos  intentando  resolver  es  homogénea  de  grado  a  respecto 
de  y, y1, . . .  ,y(n). 

En  primer  lugar,  es  claro  que,  cuando  a  >  0,  la  función  constante  y  —  0  (para 
la  cual  y'  —  y"  —  ■■■  =  0)  es  solución  ya  que,  efectivamente,  F(x,  y,  y', . . . ,  ?/n))  — 
F[x,  0,  0, . . . ,  0)  =  0 aF(x,  1, 1, . . . ,  1)  =  0. 

Para  hallar  las  demás  soluciones,  tomemos  una  nueva  función  z  dada  por  y'  =  yz. 
Obviamente,  esto  es  equivalente  a  decir  y  —  exp (f  zdx),  puesto  que 


dy 

y 


—  zdx 


log  y  = 


z  dx 


En  cualquier  caso,  si  derivamos  sucesivamente  la  expresión  y'  —  yz,  obtenemos  y"  — 
y'z  +  yz'  —  yz2  +  yz'  =  y(z2  +  z'),  y'"  —  y" z  +  2 y' z'  +  yz"  —  y(z2  +  z')z  +  2 yzz'  +  yz"  = 
y(z2  +  3 zz'  +  z"),  etcétera.  Es  importante  destacar  el  hecho  de  que  siempre  aparece  una 
relación  del  tipo  y^k>  —  ygk{z,  z' , . . . ,  z^k~x^).  Entonces,  sin  más  que  sustituir,  la  E.  D. 
queda 

F{x,  y,  yz,  y(z2  +  z'),...,  ygn(z,  z', . . . ,  z(n_1)))  =  0, 
de  donde,  extrayendo  A  =  y, 

yaF(x,  1,  z,z2  +  z',...,  gn(z,  z' z(n_1)))  =  0, 


y  por  tanto,  como  estamos  suponiendo  que  y  no  es  la  función  nula,  el  segundo  factor  habrá 
de  ser  cero.  Claramente,  esto  puede  ponerse  en  la  forma 

G(x,  z,  z' , . . . ,  z(-n_1-))  =  0, 

que  es  una  ecuación  de  orden  n  —  1.  Si  logramos  resolverla,  tendremos  que  sus  soluciones 
serán  z  —  (¡>{x,  C\, . . . ,  (7n_i),  una  familia  dependiente  de  n  —  1  constantes.  Entonces,  las 
soluciones  de  la  E.  D.  original  serán 


y  =  exp  (J  <!>{x,Ci, . .  .,Cn- 1)  dxj  , 
lo  que,  al  integrar,  introduce  la  n-ésima  constante.  Así,  puede  ponerse 


y  =  Cn  exp 


4>{x,  Ci, . . . ,  Cn-i)  dx 


RECETA  13.  Si  la  ecuación 

F{x,y,y',-  •  •  ,y(n 

--  0 

es  tal  que,  para  a  fijo,  F  cumple 

F(x,  Xuo,  Xu\, . . . ,  Xun 

)  =  XaF(x,  uq,  u\ 

,... 

7  un)t 

entonces  el  cambio  de  función 

dado  por  y'  = 

yz 

(es  decir,  y  — 

exp (f  zdx))  hace  que  el  orden  se 

reduzca  en  uno. 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


Ejemplo  13.  Resolver 

3  x\{y')2-yy")  =  y2. 

Esta  ecuación  es  homogénea  de  grado  2  respecto  de  y ,  y' ,  y"  pues,  al  sustituir  y,  y' ,  y" 
por  A  y,  Xy' ,  \y" ,  la  ecuación  queda  multiplicada  por  A2.  Para  resolverla,  hacemos  el  cambio 
y  por  ^  dado  por  y  =  exp (f  z(x)  dx ),  con  lo  cual 


✓  =  (/_-*),  y»=(^)exp  (/,*). 


Sustituyendo  en  la  E.  D.  queda 


3aú 


z2exp  2  /  zcte  )  —  (V  +  z2)exp  2  ¡  zdx 


=  exp  (2  /  z  dx  )  , 


de  donde,  simplificando  exp(2  J  2  dx),  se  sigue  —3 x2z'  =  1,  ecuación  en  variables  separadas. 
Tras  ponerla  como  dz  =  —^x~2  dx  la  integramos,  con  lo  que  z  =  |aT  1  +  C\ .  Por  último, 


y  —  exp  í  /  z(x)  dx  \  —  exp  (  /  ( 


<3X 


exp  Q  log  x  +  C\x  +  C‘¿)  =  x1/3 


1  +  Ci)  dxj 
exp  (Cía;  +  C2), 


luego  la  solución  de  la  ecuación  original  es,  sin  más  que  cambiar  la  notación  de  las 
constantes,  la  familia  de  funciones  y  =  K\x1^eK’lX . 


Pequeña  bibliografía  en  castellano 


F .  Ayres,  Ecuaciones  diferenciales,  Col.  Schaum,  McGraw-Hill,  México,  1988. 

M.  Braun,  Ecuaciones  diferenciales  y  sus  aplicaciones,  Grupo  Editorial  Iberoamérica, 
México,  1990. 

M.  de  Guzmán,  Ecuaciones  diferenciales  ordinarias.  Teoría  de  estabilidad  y  control, 
Alhambra,  Madrid,  1975. 

A.  Kiseliov,  M.  KrasnovyG.  Makarenko,  Problemas  de  ecuaciones  diferenciales  ordi¬ 
narias,  3.a  ed.,  Mir,  Moscú,  1970. 

F.  Marcellán,  L.  CasasúsyA.  Zarzo,  Ecuaciones  diferenciales.  Problemas  lineales  y 
aplicaciones,  McGraw-Hill,  Madrid,  1990. 

R.  K.  NagleyE.  B.  Saff,  Fundamentos  de  ecuaciones  diferenciales,  2.a  ed.,  Addison- 
Wesley  Iberoamericana,  Wilmington  (Delaware,  USA),  1992. 

S.  Novo,  R.  Obaya  y  J.  Rojo,  Ecuaciones  y  sistemas  diferenciales,  AC,  Madrid,  1992. 

G.  F.  Simmons,  Ecuaciones  diferenciales  con  aplicaciones  y  notas  históricas,  2.a  ed., 
McGraw-Hill,  Madrid,  1993. 

D.  E.  Zill,  Ecuaciones  diferenciales  con  aplicaciones,  2.a  ed.,  Grupo  Editorial  Ibero¬ 
américa,  México,  1988. 


47 


APÉNDICE:  Métodos  básicos 
para  calcular  integrales  indefinidas 


l.°  Reglas  principales  de  integración: 

(1)  Si  F'(x)  =  f(x),  entonces 

I  f(%)  dx  =  F(x)  +  C, 
con  C  una  constante  arbitraria. 

(2)  J  af(x)  dx  —  a  j  f(x)  dx.  donde  a  es  una  constante. 

(3)  J (f(x)±g(x))dx  =  J  f(x)dx±  J  g(x)dx. 

(4)  Si  j  f(x)  dx  —  F(x)  +  C  y  u  —  <p(x),  entonces 

í  f(u )  du  —  F(u)  +  (7, 


siendo  du  =  <p'(x)  dx.  En  particular,  si  F'(x)  =  f(x),  se  tiene 

í  f(ax  +  b)  dx  —  -F(ax  +  b)  +  C,  a  ^  0. 

/  a 


2.°  Tabla  de  integrales  inmediatas: 

/q.k+1 

xk  dx  —  - - f  C,  k  —1. 

k  +  1 


f  clx 

(b)  /  —  =log|x|+C. 

J  X 


(c)  /  ex  dx  =  ex  +  C. 
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Métodos  clásicos  de  resolución  de  E.  D.  O. 


(d)  /  ax  dx  — - b  C,  a  >  0. 

J  log  a 

(e)  í  sen  x  dx  —  —  eos  x  +  C. 


m  J  eos  x  dx  —  sen  x  +  C. 
í  dx  í , .  o 


eos2  x 


—  /  ( 1  +  tg2  x)  dx  —  tg  x  +  C . 


sen2  x 


—  /  (1  +  cotg2  x)  dx  =  —  cotg  x  +  C. 


dx  „  /X  7T\ 

- =  log  seca;  +  tga;  +  C  =  log  tg  -  +  -  +  C. 

eos  x  V  2  4  / 

dx  |  ,  a; 

- =  log  cosec  a;  —  cotg  a;  +  G  =  log  tg  —  +  G . 

sena;  2 


(k)  sh  x  da;  =  ch  x  +  C. 

(l)  í  ch  x  dx  —  sh  x  +  C. 


dx  _ 

— 2 —  =  thx  +  G. 
clr  x 


dx 

sh2  x 


2 —  =  —  coth  a;  +  C. 


dx  1  a;  1  a; 

; - 77  =  -  arctg  — b  G  = - are  cotg  — b  Gi,  a  ^  0. 

2  +  a2  a  a  a  a 


da;  1  ,  a;  —  a 

~2 - 2  =  h"~  loS  ~~T~  +  C,  a  0. 

ar  —  2a  x  +  a 


dx  1  a  +  x  1  ,  a;  .  ^ 

- 5-  =  —  log  -  +  C  —  -  arg  tli  -  +  G,  a  /  0. 

az  —  xz  2a  a  —  x  a  a 


/x2  +  a 


!x2  —  a 


=  log  a;  + 


/a2  —  ad 


■  arg  sh  —=  +  C. 
Va 

a  >  0 

arg  ch  —=  +  C. 
Va 

a  >  0 

—  -b  C i ,  a  >  0. 

a 

Integrales  indefinidas 
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3. °  Método  de  sustitución: 

Si  en  la  integral  j  f(x )  dx  se  efectúa  la  sustitución  (o  cambio  de  variable )  x  =  <p(t ) 
se  tiene 

J  f(x)  dx  —  J  f((p(t))(p'(t)dt. 

La  función  cp  se  procura  elegir  de  tal  manera  que  el  segundo  miembro  de  la  fórmula  anterior 
tenga  una  forma  más  adecuada  para  la  integración. 

4. °  Integración  por  partes: 

Si  tenemos  dos  funciones  u  =  p(x)  y  v  —  ip(x),  se  verifica 


/ 


udv  =  uv  — 


I 


v  du , 


donde  du  —  p'(x)dx  y  dv  —  ip'{x)dx.  En  este  método,  u  y  v  hay  que  buscarlas  de  tal 
forma  que  la  nueva  integral  que  aparece  sea  más  asequible  que  la  de  partida. 


